
X-645-70-4§8 

preprint 


A GENERAL COMPUTER DATA 
PROCESSING SYSTEM: DOCUMENTATION 

STATION 

MAGNETOMETER PROGRAM 


H. J. GILLIS 

r ^ ' » . ' - I v. . 


H0NV89 JLRjdNIi 
.lllOVi I1S VSVM 
Q3A!3D3y 


lasi aw 


NOVEMBER 1970 


GODDARD SPACE FLIGHT CENTER 

' GREENBELT, MARYLAND 


82 


(CODE] 


(CATEGORY) 


X-645-70-458 


A GENERAL COMPUTER DATA PROCESSING SYSTEM: 
DOCUMENTATION OF THE ATS-5 GROUND STATION 
MAGNETOMETER PROGRAM 


H. J. Gillis 


November 1970 


NASA GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 20771 




CONTENTS 


Page 

l. INTRODUCTION 1 

1. General Data Operations in this Program 1 

2. Functions of this Documentation 3 

H. COMPUTER INFORMATION 4 

m. DATA DESCRIPTION AND SOURCE 5 

IV. PROGRAM PURPOSE AND OUTPUTS 5 

V. GROUND STATION DATA TAPE FORMAT . . 6 

VI. SPECIAL SITUATIONS ON DATA TAPE 7 

VII. STEPS TAKEN TO ACCOMMODATE SPECIAL SITUATIONS 

ON GROUND STATION DATA TAPES 8 

VRI. STATION, ADDITIONAL EXPERIMENT CHANNEL (R), 

AND DATA YEAR ID CODES 10 

IX. DATA TIME VALUE ASSIGNMENT 11 

X. MAGNETIC FIELD MONITOR DESIGN DATA SAMPLE 

TIME INTERVALS 12 

XI. PLOT B (NON- AVERAGED DATA DISPLAY) FORMAT ..... 12 

XII. PLOT C (AVERAGED DATA DISPLAY) FORMAT 13 

XIII. EQUATIONS USED IN PROGRAM. 14 

XIV. SPECIFICS OF PROGRAM LOGIC 14 

1. Main Program 15 

2. Subroutine MSZRDP 16 

3. Subroutine MSC LDP 17 

4. Subroutine ATSGRT. 17 

5. Subroutine PICK. 18 


iii 



CONTENTS (Continued) 


Page 

6. Subroutine ATSGPB 19 

7. Subroutine ATSGPR 21 

8. Subroutine AT SGPC 21 

9. Subroutine ERRSET 23 

10. Subroutine PLTND 23 

11. Subroutine EXIT 23 

12. Program Input/ Output Details 23 

A. Inputs 24 

1) the data tape from the ATS-5 ground station 24 

2) the TI, TF selector program input cards 24 

B. Outputs 26 

XV. ACKNOWLEDGEMENTS 27 


iv 



APPENDICES 


Appendix Page 

A CONCEPT-LEVEL PROGRAM FLOW CHART A-l 

B FORMAT OF THE ATS-5 GROUND STATION 

MAGNETOMETER DATA TAPE B-l 

C PROGRAM LISTING C-l 

D "AUTO-FLOW" PROGRAM FLOW CHART D-l 

E SAMPLE OF THE NUMERICAL DATA VALUE 

PRINT-OUT DISPLAY OUTPUT E-l 

F SAMPLE OF THE MICROFILM NON-AVERAGED 

DATA DISPLAY OUTPUT (PLOT B) F-l 

G SAMPLE OF THE MICROFILM AVERAGED 

DATA DISPLAY OUTPUT (PLOT C) G-l 

H DEFINITION OF THE H, D, Z DATA 

COORDINATE SYSTEM H-l 

I IBM 1800 PROGRAM FOR COPYING ATS-5 

GROUND STATION TAPES 1-1 

J LISTING SHOWING THE IBM 360 JOB CONTROL 

CARDS AND THE TI, TF CARDS FOR RUNNING 
THE ATS-5 GROUND STATION MAGNETOMETER 
DATA PROCESSING PROGRAM J-l 


v 



A GENERAL COMPUTER DATA PROCESSING SYSTEM: 
DOCUMENTATION OF THE ATS-5 GROUND STATION 
MAGNETOMETER PROGRAM 


INTRODUCTION 

A computer program capable of transforming and displaying specific time- 
dependent data recorded on an indefinitely large number of tape reels is a 
common need within scientific activities. The ATS-5 Ground Station Magnetome- 
ter Data Processing Program described in this documentation is of this type. It 
is suggested that the data operations accomplished by this processing system 
are of sufficiently general application that the system could be used to process 
other types of data with little modification. To enable a determination of the 
applicability of this system for a particular data processing use, a list of these 
data operations is presented below. Also, a concept-level "plain language" flow 
chart for their implementation is presented in Appendix A, and a commented 
listing of the production program is presented in Appendix C . The program 
logic is summarized under "Main Program" on page 15, and in the flow chart in 
Appendix A. This documentation uses an original instructional device involving 
numbered comments to aid the reader to understand the program. This is ex- 
plained more fully on pages 4 and 14. 

The effectiveness with which the program has processed the ATS-5 data suggests 
that it successfully avoids the weak points of too much or too little sophistication 
which has troubled other data processing programs observed by the author. Data 
systems written by the author using the techniques described in this publication 
have kept the number of unprocessed ATS-5 flight and ground data tapes on hand 
to zero and are concise enough to be completed in four months by one experienced 
person. The production program card deck will be supplied upon request. 

1. General Data Operations in this Program 

A. Expression of data time in terms of a single physical unit (milliseconds) 
for computational convenience for data extending over minutes, hours, 
days, or years (see explanation under "Data Time Value Assignment" 
on page 11). 

B. Selection of time intervals of data to be processed by designation of 
their beginning times (TI) and end times (TF) on program input cards 
in chronological order. 

C . Search through the tape to data with times within the time interval cur- 
rently selected for processing. 
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D. Confirmation that the tape used in the computer run is the one in the 

data tape library bearing the data time intervals selected for processing 
by the TI, TF program input cards. This is done by print-out of the 
beginning and end time of the data tape mounted on the tape drive which 
shows if this tape contained the data time interval currently selected 
for processing as indicated by print-out of its TI and TF. This program 
proceeds to the next TI, TF selector card if the mounted data tape post- 
dates the first TI, TF selector card. It terminates the computer run if 
the end of the data tape is encountered while searching for or processing 
data in the currently selected time interval. 

NOTE: In other programs using this general data processing scheme such as 

those for ATS-5 and OGO-6 satellite flight data which, contrary to this 
program, are intended to process many data tapes in a single computer 
run differing action is taken when the mounted tape does not contain 
data in the currently selected time interval. Instead, the computer 
operator is provided with on-line printer messages showing the begin 
and end time of each tape he mounts enabling him to ultimately mount 
the tape containing data in the currently selected time interval (also 
printed on-line) if this tape has been made available to him, and through 
some oversight he at first mounts some other tape. These programs 
permit the operator to proceed to the next TI, TF selector card if the 
proper tape is not available. Also, programs using this general data 
processing scheme avoid the considerable programming complication 
of a "data time span versus relevant data tape identification number" 
directory logic module for automatic on-line printer tape mounting re- 
quests. In these programs this directory function is achieved by desig- 
nation of the data tape identification numbers in the sequence in which 
they are to be used in the run by the person submitting the program. 

This person refers to the data tape shipping slips for tape begin and 
end time information to do this. This scheme has been found practicable 
even when as many as 12 or more ATS-5 satellite PB data tapes have 
been submitted in a computer run directly producing output data dis- 
plays. The sheet of on-line printer messages showing begin and end 
times of the time intervals selected for processing, the begin and end 
times of the data tapes actually mounted, and printed indications of the 
computer operator's response to the respective tape mount requests 
constitutes a report that would show any deviations from the data tape 
utilization instructions submitted with the job. The only prerequisite 
on the data organization on a tape for this general data processing 
scheme to process all data on it is that the data within any given file be 
in chronological order — the various files of the tape do not have to be. 
Also, it is permissible for data missing from a tape to appear on another 
tape. 
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E. Ignore any data on the tape with time which is out of chronological 
order. 

F. Do processing unaffected by any missing data, i.e. irregular jumps in 
data times. 

G. Unpack the data values from their fields of various bit lengths in each 
logical record (a grouping of bits on tape comprising the data in a 
definite format which by repetition forms the data tape) into fortran 
variables for processing in the main program. 

H. Mathematical conversion of the data into appropriate physical units. 

I. Isolation of collections of data {referred to as Plot B arrays) of specific 
time intervals relative to the beginning of the day of data to do statistics 
on and/or to graphically display on microfilm as a function of time 
(termed the Plot B display) . 

J. Obtain collections of time- averaged data points (referred to as Plot C 
arrays) to produce "quick-look*' summarizing graphical microfilm plots 
of events over specific time intervals relative to the beg innin g of the 
day of data (termed the Plot C display). These "summary" microfilm 
frames appear right after the end of the microfilm frames which they 
summarize. The time interval over which an average is computed is 

a variable of the program presently set at 30 seconds. 

K. In graphical microfilm display, choice of ordinate and abscissa scale 
(independently) to display data plotted in any given frame with the best 
resolution. The choice of scale is based on the prevailing data sampling 
time interval for the abscissa scale, and the range in data values plotted 
in the microfilm frame for the ordinate scale. 

L. Display numerical values of processed data as a function of time on 
paper print-out. 

2. Functions of this Documentation 


A. Define the type of data processing accomplished by this program pointing 
out the methods used, and their usefulness for writing programs to do 
similar processing on data of other types (see Appendix A for flow chart 
illustrating how various general data operations are done in this program.) 
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B. Describe the data, source of data, and data tape format specifically 
processed by this program (see Appendix B for the data tape and logical 
record format) , 

C. Itemize all subroutines and their functions. 

D. Explain in detail the program logic using numerous plain language 
comments prefacing the program statements to which they apply (see 
program listings in Appendix C and detailed "Autoflow" flow chart in 
Appendix D) . These comments are numbered to identify corresponding 
specific program statements in the program listing by comment number 
reference in explanations of program logic in the written documentation 
and in the flow charts in Appendices A and D. 

E. Describe the various display outputs of the program and show samples 
of them (see Appendices E, F, and G). 


COMPUTER INFORMATION 

The program is written for the IBM/360, Operating System release number 16 
or higher with the new SCORS Stromberg Carlson 4020 plotter subroutine pack- 
age. However, the main program can be used without change for any plotter since 
all plotting logic is "modular", i.e., done by subroutines called by the main pro- 
gram. To use the program to drive a plotter other than the SC4020, it is neces- 
sary only to replace the SC4020 subroutine calls by calls to the subroutine pack- 
age of the replacing plotter. 

The main program and all but one of the subroutines are in Fortran IV. A small 
subroutine named PICK used for unpacking the data values from the bit sequences 
in the logical records on the data tape is in IBM/360 Assembly Language. Sub- 
routine PICK is sufficiently general to be used for data "unpacking" in a logical 
record of any format. 

All program input and output operations are done by Fortran IV in the main pro- 
gram since IBM/360 Fortran has sufficiently versatile capabilities to read data 
tapes of arbitrary format without the need for the less simple I/O instructions 
of the assembly language. 

The amount of IBM/360 core required by the program is 200,000 bytes. The 
IBM/360 computer time required on the Model 75 is about 20 minutes for each 
day of data. On the IBM/360 Model 91 the program runs twice as fast. 
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DATA DESCRIPTION AND SOURCE 


The data read from a tape and processed by this program is raw magnetometer 
readings of the components of the geomagnetic field in the H, D, Z coordinate 
system (see Appendix H for definition) and readings from an additional experi- 
ment channel (called R readings) at a specific location on the earth's surface as 
a function of time. This location is at a Canadian geophysical ground station 
situated near the foot of the geomagnetic line of force passing through the ATS-5 
satellite. Two such stations are operating: LYNN LAKE and THOMPSON in 
Manitoba, Canada. All of the data on a given tape will be from one of these sta- 
tions. The output displays of the program are identified by labelling showing 
the name of the station to which they apply. 

At both stations, data is written on the tapes by inexpensive tape decks fed by 
the magnetometer. The tape decks are specially designed to directly write the 
data on tape in physical records (a division of data bits on tape containing special 
non-data control and check bits and bounded on both sides by "inter-record gap" 
tape marks) in accordance with IBM/360 tape physical record specifications. 

This enables reading of the data tapes by the IBM/ 360 without intermediate tape 
format conversion procedures. Approximately 30 days of data is written on a 
tape, the tape deck-magnetometer digital data recording unit being unattended 
for this length of time. 

The phrase "set of H,D, Z,R data" used throughout this documentation signifies 
the values of the H,D, Z, and R data readings- sampled at the same instant of 
time and recorded on tape in the normal course of operation of the ground sta- 
tion digital magnetic data recording unit. 


PROGRAM PURPOSE AND OUTPUTS 

The purpose of the program is to convert the raw H,D, Z coordinate system 
geomagnetic field component data count readings and auxiliary experiment chan- 
nel R readings (if channel R is being used for magnetic data) on tapes from 
ATS-5 ground stations into gamma (equal to 10 -5 gauss) units for a selected 
time interval of interest, and to display the results. 

The display output of the program is as follows: 

1. time-ordered numerical print-out of H,D, Z (in gamma units), and the 
name of the experiment on the R channel and the R data values (see 
sample in Appendix E). 
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2. graphical plot of the instantaneous H,D, and Z values plotted individually 
as a function of time on Stromberg Carlson 4020 plotter microfilm 
(called Plot B, see sample in Appendix F). 

3. graphical plot of 30-second average H, D, and Z values (in gammas) 
plotted individually as a function of time on SC 4020 plotter microfilm 
(called Plot C, see sample in Appendix G), 

NOTE: The time assigned to each averaged H, D, or Z value is the mid-time of 

the 30 -second time interval containing the data used for calculating the 
respective average. 


GROUND STATION DATA TAPE FORMAT 


The designed tape format is as follows (see Appendix B for detailed schematic 
diagram of the tape format) : 


1. Physical record (alternatively called "block") control bits and inter- 
record gap length — within IBM specifications. 

2. Physical record length — 7200 IBM/360 bytes. 

3. Logical record characteristics (see Appendix B for detailed description 
of a logical record). 

A. Length - 72 IBM/360 bytes 

B. Logical records per block — 100 

C. Data contents of a logical record (in order of occurrence on tape) 


1) time (day of year, hour, minute, second) of first of the ten sets 
of H,D, Z,R data counts in the logical record 

2) station ID code 

3) R experiment channel code 

4) data year code 

5) first of the ten sets of H, D, Z, R data counts in the logical 
record 
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6) 48 dummy bits not presently used 

7) the nine remaining sets of H,D, Z,R data counts of the logical 
record 

D. Time span covered by logical record — depends on the time in- 
terval between sampling of each set of H,D, Z,R data. This time 
interval is set manually by a control on the Magnetometer Field 
Monitor (see "Design Sample Time Intervals" below). The data 
sample time interval is usually one second, 

E. Data files per tape — 1 
SPECIAL SITUATIONS ON DATA TAPE 

The actual tapes differed from the original design specifications in some blocks 
as follows: 

1. block length less than 7200 bytes (some as short as one byte) 

2. control bits not within IBM specifications 

3. data sampling time interval different from design value 

4. year code (see "ID Codes" below) for 1969 wrong (for entirety of tape) 

5. day of year of a logical record occasionally wrong 

6. occasional blocks with incorrect format in that the block does not begin 
with the beginning of a logical record (the data time field) causing all 
following logical records of the block to also be off-format 

7. spurious extra ends of file on data tape 

The IBM/360 reacts to short blocks (less than 18 bytes) and/or incorrect con- 
trol bits as though an I/O error for the block involved occurs. This causes 
multiple re-reads of the block by the IBM/360 rapidly wearing through the data 
tape. 
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STEPS TAKEN TO ACCOMMODATE SPECIAL SITUATIONS 
ON GROUND STATION DATA TAPES 


1. A computer that re-reads blocks with faulty control bits and/or length 
less than 18 bytes zero times was found. A program for it, shown in 
Appendix I, was written to copy blocks of the original tape padded out 
to 7200 bytes by hex 9's if short and with a hex 1 inserted in the last 
hex digit of the dummy bit field of the first logical record as an indi- 
cator if read with I/O error. All data of a block with such an indicator 
is flagged with an "F" in the microfilm display and numerical print-out 
since data read with an I/O error may be unreliable. The copied data 
tape has no faulty control bit or short block problem when run on the 
IBM/360. This computer used for copying is the IBM/1800 located in 
Goddard Space Flight Center, Building 2 ground floor. It can be oper- 
ated by the programmer himself. 

2. The program calculates the data sampling time interval for the data of 
a logical record (with valid time information) by dividing by 10 the 
difference between the first data sample time (the first 9 hex digits of 
the logical record) and that of the next consecutive logical record if it 
has valid time information. This handles any operationally imposed 
data sampling time interval changes within the block. Tests deter- 
mining what constitutes valid times are described in (3.) following. 

The 10 sets of H,D, Z,R data of any logical record having a data 
sampling time interval differing from design specifications are pro- 
cessed as usual, but are flagged with a "T" in the microfilm display 
and the numerical print-out. 

3. A logical record which has an invalid time field (first 9 hex digits) or 
which is followed by a logical record having an invalid time field 
(making determination of data sampling time interval impossible) is 
not processed. Following are conditions defining an invalid time field 
of a logical record: 

A. day of year, hour, minute, or second having an impossible value, 
e.g., day of year greater than 366 or less than 1 (if this occurs the 
time field of the logical record is set to zero as an indicator in- 
stead of the equivalent number of milliseconds since the "ZERO 
YEAR" time origin as described below). 

B. time of the logical record greater than time of following logical 
record (chronological order test). 
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C. time difference between logical record and the following one 
greater than maximum value (10 seconds) prescribed by design 
specifications (guards against the use of an incorrect data sampling 
time interval in assigning a time to each of the 10 sets of H,D, Z,R 
data in the logical record) . 

4. To enable processing of the data in the occasional short blocks (less 
than 7200 bytes) on the original tape the program puts hex 9's in that 
part of the fortran array (IDAT) receiving the data block which is not 
filled by the block due to its shortness. In this way data from a previ- 
ous block will not remain at the end of the array to be incorrectly in- 
terpreted as data of a new block. Thus, the program will ignore any 
blocks having hex 9’s in the time field of the 2nd logical record or 
before. Also, any logical records having hex 9's in it or in the time 
field of the following logical record is recognized as the last logical 
record of a short data block and is not processed. 

5. In the absence of the correct year code field in all logical records of 
data tapes of year 1969, (year code for '69 as well as '70 was zero), 
the change of day of year to a value less than 365 was used to indicate 
beginning of 1970 data. Due to the unreliability of the year code field 
in the data tapes, the correct year at the beginning of the tape being 
processed is supplied by means of a fortran statement placed in the 
main program. 

6. Data tapes with spurious ends-of-file are processed by means of a 
specially modified program deck which considers the tape as one having 
multiple IBM/360 data sets and uses the required special IBM 360 job 
control cards. The actual number of files on the tape must be entered 
into this special program which is not shown in this publication. A copy 
of this program can be obtained from the author. 

7. Avoidance of data loss in tape blocks with format errors was accom- 
plished by use of special program logic. This logic searches for the 
first occurrence in the tape block of a distinctive bit configuration (the 
dummy bit field of the logical record as shown in Appendix B). Since 
this relatively easy to locate bit configuration occurs at a fixed position 
within a logical record, the hexadecimal digit beginning the first logical 
record in the tape block is also determined once the first dummy bit 
configuration is found. If the hexadecimal digit beginning the first 
logical record is not the first hexadecimal digit in the tape block, the 
fortran array containing the data of the tape block is shifted so that it 
is before processing continues. This format error corrector fails if 
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the logical record length is not constant as occasionally happens. The 
program utilizing this corrector logic is not shown in this publication. 
Copies of this program can be supplied however. 


STATION, ADDITIONAL EXPERIMENT CHANNEL (R), 

AND DATA YEAR ID CODES 

See Appendix B for where these codes appear in each logical record of the data 
tape. The field length for each code is 4 bits, i.e., one hexadecimal digit. The 
codes, in hexadecimal, follow: 

1. Station Codes 

A. 1 = LYNN LAKE 

B. 2 = THOMPSON 

C. 3 = WINNIPEG 

D. 4 = THE PAS 

2. Additional Experiment Channel (R) Codes 

A. 0 = not used 

B. 1 = H magnetic field component 

C. 2 = D magnetic field component 

D. 3 = Z magnetic field component 

E. 4 = proton experiment (total magnetic field) 

F. 5 = other experiment (photometer, etc.) 

3. Year Codes 

A. 9 = 1969 

B. 0 = 1970 

C. 1 = 1971 
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D. 2 = 1972 


E. 3 = 1973 
etc. 


DATA TIME VALUE ASSIGNMENT 

Each logical record has 10 sets of H, D, Z, R data each set corresponding to an 
instant of time. Only the time (in day of year, hour, minute, second) of the first 
set of H, D, Z, R data is given in the logical record. The time for each of the 
9 remaining sets of H, D, Z, R data is obtained by adding to the time of the first 
H, D, Z, R set the data sampling time interval (between each H, D, Z, R set of 
data) the number of times appropriate for the particular set of H, D, Z, R to 
which a time value is being assigned. For example, to get the time of the third 
H, D, Z, R set the data sampling time interval is added to the time of the first 
H, D, Z, R set twice. As mentioned above, the data sampling time interval for 
the logical record is the difference between the time of its first set of H,D, Z,R 
data and the time of the first H, D, Z, R set in the next consecutive logical 
record divided by 10. If either of these two times do not pass the time validity 
tests described above, no data sampling time interval is calculable and the data 
of the logical record is ignored. This method enables automatic handling of 
change of the operating sampling time interval wherever it occurs on the data 
tape. 

To do computations with time, such as time addition to get data point times, or 
taking time difference between begin time of a microfilm frame and the data 
point time being plotted (to get its abscissa in one system of units as required 
by plot subroutines), etc., any time value used in the program is converted from 
calendar units (year, day of year, hour, minute, second) to the equivalent num- 
ber of milliseconds the time value is from an arbitrary time origin. This time 
origin (called ZERO YEAR) is prior to the earliest time value in the data tapes 
so that the "number of milliseconds since ZERO YEAR" equivalent of a data 
point time will never be negative. The arbitrary point in time taken as ZERO 
YEAR is zero hour, minute and second of the first day of any year if this time 
precedes the earliest data to be processed. However, ZERO YEAR must not 
pre-date the data so much that the "number of milliseconds since ZERO YEAR" 
expressing the latest time used in the program is so large as to exceed the com- 
puter storage allotted to it. In this program, times (in "milliseconds since zero 
year") are in fortran double precision which allows (for the IBM/360) processing 
of time's near enough to ZERO YEAR so that their equivalent in "milliseconds 
since ZERO YEAR" is 17 significant digits or less. Thus many years of data 
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can be processed without change of the ZERO YEAR (set by a fortran arithmetic 
statement to 1969 at beginning of main program). 

The program uses two auxiliary subroutines for operations on time. These are 
subroutine MSZRDP to go from time in calendar units (year, day, hour, minute, 
second) to the equivalent in milliseconds since ZERO YEAR, and subroutine 
MSCLDP to go from "milliseconds since ZERO YEAR" to the more convenient 
equivalent in calendar units (year, day of year, hour, minute, second). 


MAGNETIC FIELD MONITOR DESIGN DATA SAMPLE TIME INTERVALS 

The operating data sample time interval is set by a control switch on the Mag- 
netic Field Monitor and is seldom changed. The usual data sample time interval 
is one second, i.e., the time interval between times of consecutive sets of H, D, 
Z, R data is usually one second. The possible data sample time intervals by 
original monitor design are .1 sec, 1 sec, 2 sec, 5 sec, or 10 sec. 


PLOT B (NON-AVERAGED DATA DISPLAY) FORMAT 

Each microfilm frame consists of a graph of the H component (in gammas) as a 
function of time, a similar D component graph, and a similar Z component graph 
of the same frame. See Appendix F for a sample plot B. 

The time scale for the H, D, and Z component plots on the frame is the same 
and depends on the data sampling time interval for the first set of H, D, Z data 
of the microfilm frame so that the time scale optimum for the prevailing sam- 
pling time interval at the start of the plots is used. The time scale (the data 
time length covered in the microfilm frame) that we considered optimum as a 
function of data sampling time interval is as follows: 


DATA SAMPLING 

DATA TIME LENGTH 

TIME INTERVAL 

COVERED IN FRAME 

.1 second 

1 minute 

1 second 

6 minute 

2 second 

12 minute 

5 second 

1 hour 

10 second 

1 hour 
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The program causes the begin time and end time of any plot B to delimit an 
integral multiple of the chosen data time length from time zero of the day of 
data, and to include at least one data point positioned properly on the time scale 
even though data gaps may exist. This facilitates data comparisons for different 
days. A vertical line appears on the microfilm frame every sixth of the se- 
lected data time length covered by the frame. The time intervals between these 
vertical lines are further sub-divided by 10 short vertical "tic" marks. 

The vertical scale for the H (in gamma units) plot of the frame depends on the 
maximum and minimum values of the H component values being displayed in the 
plot. Similarly for the vertical scale of the D plot, also the vertical scale of 
the Z plot. Only six vertical scales, i.e., low and high limit of vertical scale 
for a given component are possible. The specific vertical scale chosen for the 
component plot on the particular microfilm frame is the first of the following 
six scales which includes both the maximum and the minimum of the component 
values displayed in the plot. 

VERTICAL PLOT SCALES 

(the one displaying the 
component with the best 
resolution is chosen for 
the microfilm frame) 

-60. to +60. 

-150. to +150. 

-300. to +300. 

-600. to +600. 

-1200. to +1200. 

-2400. to +2400. 

A horizontal line appears on the respective component plot on the microfilm 
frame every sixth of the data range (in gamma units) covered by the selected 
vertical scale. The scale between these horizontal lines is further sub-divided 
by 10 short horizontal "tic" marks. 


PLOT C (AVERAGED DATA DISPLAY) FORMAT 

A Plot C appears every time data spanning an hour of time has been displayed 
by plot B's. Its purpose is to summarize (as 30-second time interval averages) 
for each individual component the instantaneous data displayed in the intervening 
(an "hour's worth") plot B's. This provides a "quick-look" capability to the 
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scientist using the microfilm. For example, a plot C summarizes and appears 
after every ten plot B's when the plot B data time length is six minutes (which 
is the norm) . 

The scheme for the vertical and horizontal scales for plot C are the same as 
those for plot B described above except the horizontal scale (the time scale) is 
non-variable and always covers one hour of data. See Appendix G for a sample 
of plot C. 


EQUATIONS USED IN PROGRAM 

1. To obtain the H, D, or Z field components values from magnetometer 
"counts" for the individual component as recorded on the data tape. 

FIELD COMPONENT (in gammas) = COUNTS FOR COMPONENT x .976408 - 2000. 

2. To provide reliability that a 30-second average component value does 
approximate the true field component value at the time assigned it, no 
average component value is established for the 30-second interval 
unless at least 1/3 of the maximum number (30 seconds divided by the 
data sampling time interval) is available in the 30-second interval for 
computing its average. Thus, the minimum number of component values 
required (CPCTMN) for getting a 30.- second average component value is. 
given by 


_ J. 30000. _ 10000. 

CPCTMN 3 • TSpLST TSPLST 

where TSPLST is the data sampling time interval at end of the 30- 
second interval. TSPLST rarely changes. 


SPECIFICS OF PROGRAM LOGIC 

Only the general features of the program are discussed here since in-depth 
details intended to expedite the programmer’s learning task appear as numbered 
comments in the main program listing in Appendix C. Numbering of the com- 
ments makes it possible to locate the specific group of Fortran statements being 
explained in the text below which cites comment numbers, or being represented 
by flow chart blocks (Appendix A has general flow chart. Appendix D has detailed 
"Autoflow" flow chart) adjacent to which the corresponding listing comment 
statement numbers appear. 
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1. Main Program 


The scheme of data processing used in the main program resulted from 
attempting to attain a specific goal: to write the most brief, straightforward, 
efficient program for processing the given data. The ability to process only 
selected time intervals (by designating their begin and end times, i.e., TI, TF's 
on data cards read by the program) or the entire data tape was also included in 
the program. The general data processing scheme used here for ATS-5 ground 
station data has also been used quite successfully with large volumes of data 
from other sources (ATS-5 PCM and PFM flight data, and also OGO-6 flight data). 

The data processing steps done by the main program are: 

A. read next TI, TF or end-of-computer run indicator (XX) on data card 
(comment 2); the XX card is always last data card. 

B. read next block on data tape until last time of block is greater or equal 
to TI (comments 2 to 10); terminate run if end-of-tape encountered in 
this "search mode". 

C. get begin time of block accessed by previous step; back to step "A" if 
this begin time is greater than TF, i.e., beyond present TI, TF (com- 
ment 10). 

D. process the logical records of the block that are within the present TI, 

TF by storing the time and the H, D, Z (in gammas) and R of each of 
the 10 sets of H, D, Z, R counts of each logical record in the arrays 
("B arrays") containing the data to be displayed by a plot B (comments 
11 to 34). 

NOTE: The extra step of storing data in an array instead of plotting it directly 

is necessary in order to pre-determine the best vertical scale for the 
H, D, and Z plots individually before doing the plot B or the plot C. 

During this processing, 3 situations may occur: 

1) data time larger than TF encountered — in this event go to step "A" 
(comment 16). 

2) end of block on data tape reached in processing — in this event to 
to step "B" (comment 15). 

3) data time larger than end time of present plot B encountered — in 
this event go to step "E", etc., below (comment 22). 


15 



E. get 30 second average data points for plot C from the B array before it 
is filled by data for the next plot B; store these averaged data points in 
an array termed the "C array” (comments 22 to 31). 

F. if, in step "E", a 30 second average data point with time greater than 
end time of present plot C is encountered, use the data in present C 
array to do a plot C. After this plot C is done, the C array is again 
available for storing more 30 second averaged data for the next plot C 
starting with its first location (comments 26 to 28) . 

G. now that the data in the B array has been utilized (for averaging), use 
this data to do the plot B whose end time was just exceeded; after this, 
the B array is again available. Continue storing the data from the logi- 
cal records into the B array starting with their first location (comments 
31 to 34). 

H. computer run ends by reading the last TI, TF data card which is always 
the symbol XX indicating the computer run is to be terminated (see 
step "A"), or by encountering end-of-tape (comments 2 and 7). 

2. Subroutine MSZRDP (IYR, IDY, IHR, MN, ISEC, TM) 

A. Purpose — to get "millisecond since ZERO YEAR” equivalent in double 
precision of a time known in calendar units, i.e., year, day of year, 
hour, minute, and second. 

B. Calling Sequence — 

1) IYR = last 2 digits of year (fortran integer) 

2) IDY = day of year, i.e., Julian day (fortran integer) 

3) IHR - hour of day (0 to 24, fortran integer) 

4) MN = minute (fortran integer) 

5) ISEC = second (fortran integer) 

6) TM = the equivalent in "milliseconds since ZERO YEAR" (fortran 
double precision variable) of the time specified by the 5 preceding 
calling sequence elements. 
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C. Common Section — one only, named ZROYR, to communicate value of 
ZERO YEAR, in common with main program in which value is set (for 
explanation of the term ZERO YEAR see the section above entitled 
"Data Time Value Assignment"). 

3. Subroutine MSCLDP (TM, IYR, IDY, IHR, MN, SEC) 

A. Purpose — to get the calendar units equivalent, i.e., year, day of year, 
hour, minute, second of a time known in "milliseconds since ZERO 
YEAR". Also returns month and day of month (see "common sections" 
below) . 

B. Calling Sequence — 

1) TM = the equivalent in "milliseconds since ZERO YEAR" (fortran 
double precision variable) of time specified by the 5 following 
calling sequence elements outputted by this subroutine. 

2) IYR = last 2 digits of year (fortran integer). 

3) IDY = day of year, i.e., Julian day (fortran integer). 

4) IHR = hour of day (0 to 24, fortran integer) . 

5) MN = minute (fortran integer). ■ 

6) SEC = second (fortran single precision floatingpoint variable). 

C. Common Sections — 

1) name is ZROYR, used to communicate value of ZERO YEAR to this 
subroutine, in common with main program in which value is set. 

2) name is DATE, used to return month (MNTH, typed as fortran in- 
teger but contains BCD alphabetic information) and day of month 
(IDYMTH, fortran integer) to calling program. 

4. Subroutine ATSGRT (IWD1, IWD2, IYR, TM) 

A. Purpose — to get the data time presented in the logical record in 
"milliseconds since ZERO YEAR". 
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B. Calling Sequence — 

1) IWD1 = a fortran integer word consisting of the first 4 bytes of the 
logical record, i.e., day, hour, minute and 10's digit of second (see 
format of the logical record in Appendix B) . 

2) IWD2 = a fortran integer word consisting of the second 4 bytes of 
the logical record needed only to get l's digit of seconds. 

3) IYR = last 2 digits of year of data, set by arithmetic statement at 
beginning of the calling program (the main program) instead of 
picked-up from data tape since the data year code for year 1969 is 
incorrect on the tapes. 

4) TM = the equivalent in "milliseconds since ZERO YEAR" (fortran 
double precision variable) of the data time presented in the logical 
record. 

5. Subroutine PICK (ITO, EPROM, ISW, IOFST, NRBTS) 

A. Language — IBM 360 assembly language. 

B. Purpose — to unpack a bit string of the logical record constituting one 

of its data fields and move it to a fortran integer word in the calling 

program which then uses it. 

C. Calling Sequence (all elements are integers) — 

1) ITO = address of fortran word where bits are to be moved and 
right-adjusted in this word. 

2) IFROM = address of fortran word containing the bit string com- 
prising the data field wanted. 

3) ISW = word combination switch, when non- zero allows more than 
one fortran word to be used in building the bit string being returned 
as a data field to the calling program in the fortran word specified 
by ITO above. 

4) IOFST = number of bits that contents of word specified by IFROM 
above must be shifted left in order to left-adjust the bit string to 
be moved to the calling program. 
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5) NRBTS = number of bits comprising bit string, i.e., data field 
wanted. 

6. Subroutine ATSGPB (T MB ARY, BT, ID BAR, ICT, IVSC, SCLM1, SCLM2 , 
TIPB, TFRLTH (IHSC), IF LG, ITFLG) 

A. Purpose — contains all logic of program for doing the plot B for the 
Stromberg Carlson 4020 microfilm plotter; plots all data presently 
contained in the B arrays which are communicated to it; the data in the 
B arrays is within the specific begin and end time chosen by the main 
program and suitable to the present data sampling time interval (see 
section entitled "Plot B Format" above). 

B. Calling Sequence (all input quantities to subroutine) — 

1) TMBARY = array containing time in double precision "milliseconds 
since ZERO YEAR" of each set of H, D, Z components plotted; 
dimensioned large enough (to 730) for the case in which the B array 
contains the greatest number of points to be plotted, i.e., when the 
data sampling time interval is 5 seconds with a 1 hour time length 
displayed on plot B. 

2) BT = two-dimensional array containing the individual H, D, Z 
component values to be plotted; first subscript selects component 
to plot (1 for H, 2 for D, 3 for Z); second subscript selects the H, 

D, Z set being plotted (data sampled at the same instant of time 
comprise a set). Each H, D, and Z value is in gamma units and 
fortran single precision. 

3) ID BAR = two dimensional fortran integer array containing the in- 
dividual identification codes used primarily in SUBROUTINE 
ATSGPR; its first element is used in this subroutine to identify for 
labelling purposes the ground station whose data is being plotted; 
the first subscript selects which code (1 for station, 2 for identifi- 
cation of extra experimental channel R, 3 for indication of year of 
data); second subscript selects the H, D, Z set to which the codes 
apply. 

4) ICT = actual number of H, D, Z sets ("data points") in array BT 
(see above) to plot. 

5) IVSC = array containing as one of its 3 elements the subscript of 
the vertical scale (see below) that displays the H component with 
best resolution on this plot B; similarly for D and Z. 
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6) SCLM1 = array containing the plot B lower limit for each of the 
six possible vertical scales for the individual H, D, or Z component; 
the specific lower limit for the individual H, D, or Z plot on this 
plot B is selected by the SCLM1 array subscript value stored in the 
IVSC array element corresponding to the component (H, D, or Z). 
For example, the lower limit of the H plot vertical scale is the 
SCLM1 array element with subscript equal to IVSC(l). Similarly 
for D and Z. 

7) SCLM2 = array defined similarly to array SCLM1 except SCLM2 
contains the plot B upper limits for each of the six possible vertical 
scales. 

8) TIPB = begin time in fortran double precision "milliseconds since 
ZERO YEAR" of this plot B; not necessarily a data time but it is 
the begin time of the closest integral multiple from zero instant of 
the data day of the time length selected for (see section entitled 
"Plot B Format") and displayed in this plot B that contains the first 
data time on the plot B. 

9) TFRLTH(IHSC) = the element of array TFRLTH selected by the 
value of subscript IHSC; array TFRLTH contains the millisecond 
equivalent of the various data time lengths that plot B displays as 
selected by the prevailing (at beginning of data in the corresponding 
B array) data sampling time interval. 

10) IF LG = array containing symbol for each set of H, D, Z values 
("data point") indicating if the data tape block it is in was read by 
the IBM/ 1800 copy program and/ or the actual data processing 
IBM/360 program with an I/O error. The symbol "F" is used to 
signify I/O error. Caution in interpreting such data should be used 
by the scientist. The blank symbol is used to indicate data read 
without an I/O error. 

11) ITFLG = array containing symbol for each set of H, D, Z values 
("data point") indicating if the data sampling time interval value 
for the set is equal to one of the design values (stored in array 
TINT). The symbol "T" is used to mark a data point at which the 
data sampling time interval is unexpected. Data marked by a "T" 
is not necessarily wrong. The blank symbol is used to indicate 
data at which the data sampling time interval is equal to a design 
value. 
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C. Common Section — name is DATE, used to return month (MNTH, typed 
as fortran integer, but contains BCD alphabetic information) and day of 
month (IDYMTH, fortran integer) to calling program. 

7. Subroutine ATSGPR (TMBARY, BT, R, ID BAR, IFLG, ITFLG, ICT) 

A. Purpose — this subroutine is called only if the symbol "PRT" appears 
on the TI, TF data card (see below); contains all logic of the program 

to do numerical print-out of the H, D, Z magnetic field component values 
in gamma units in chronological order for data within the currently se- 
lected TI, TF. Included in this print-out is the data of the additional 
experimental data channel (R) and the name of the experiment on chan- 
nel R. If the R data is either the H, D, or Z component it is in gamma 
units, otherwise the R data is printed out exactly as found on the data 
tape, i.e., in counts. The data printed out by a call to this subroutine 
is the contents of the present plot B array at time of the call. 

B. Calling Sequence — 

1) R = array containing all R channel data values for data times in 
present plot B array; in gamma units if R channel has magnetic 
data, otherwise in counts; in either case, fortran single precision. 

2) for explanation of other elements of calling sequence see explana- 
tions given for them under SUBROUTINE ATSGPB. 

C. Common Section — 

1) name -is DATE; communicates month (integer-type fortran variable 
MNTH containing 3-character alphabetic name of month), and day 
of month (fortran integer IDYMTH) communicated from SUBROU- 
TINE MSCLDP (see above) to this subroutine for printing out time 
information. 

8. Subroutine ATSGPC (TMAV, BTAV, ISBSTA, ICTC, IVSC, SCLM1, 

SCLM2, TIPC ) 

A. Purpose — contains all logic for doing an SC 4020 microfilm frame 
(a plot C) displaying the 30 second averaged H, D, and Z magnetic field 
component values in individual plots on the frame from data in the 
C arrays each time called. 
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B. Calling Sequence — 

1) TMAV = array containing times in fortran double precision 
"milliseconds since ZERO YEAR" of the H, D, Z 30 second average 
component values displayed in this plot C; dimensioned large enough 
(to 130) to hold maximum number of data times (at 30 second 
spacing) contained within the fixed data time length displayed by a 
plot C (1 hour). 

2) BTAV = two-dimensional array containing the H, D, and Z 30 second 
average component values to be individually displayed on plot C; 
first subscript selects component (1 for H, 2 for D, 3 for Z); second 
subscript selects the set of average H, D, Z values being individu- 
ally plotted where "set" in this usage signifies an average H, D, 

and Z value all applying to the same assigned data time (stored in 
array TMAV above). 

3) ISBSTA = ground station identifier code used to choose ground 
station name with which to label this plot C; the code used comes 
from the logical record that contains the first set of instantaneous 
H, D, Z component data values included in the B array last used 
for averaging and filling of the C arrays. 

4) ICTC = actual number of 30 second average H, D, Z component 
value "sets" in C array to be displayed on this plot C. 

5) IVSC, SCLM1, SCLM2 = see explanations for these arrays given 
under calling sequence for SUBROUTINE ATSGPB. 

6) TIPC = begin time in fortran double precision "milliseconds since 
ZERO YEAR" of this plot C; not necessarily a data time but it is 
the begin time of the closest integral multiple from zero instant of 
the data day of the fixed plot C time length (1 hour) which contains 
the time of the first 30 second average H, D, Z component values 
on the plot C. 

C. Common Section — 

1) name is DATE; explanation similar to that for common section 
DATE in subroutine ATSGPB above. 
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9. Subroutine ERRSET (IERNO, INOAL, INOMES, ITRACE, ADDUSE, IRANGE) 

A. Purpose — this is an IBM/360 Fortran Programming System subrou- 
tine. It prevents the operating system from terminating the computer 
run if multiple short blocks, i.e., of length (in bytes) less than that 
specified by the Job Control Language (see Appendix J) for the data set, 
are encountered on the ground station data tape. This "short block" 
situation does sometimes occur due to tape deck operation deviation 
from design. 

B. Calling Sequence — 

1) see explanation of this subroutine in IBM/360 System Reference 
Manual "Fortran IV (G and H) Programmer's Guide" (form number 
C28-6817-0) in Chapter entitled "Extended Error Message Facility". 

10. Subroutine PLTND 

A. Purpose — this is one of the subroutines of the programs used in 
generating the SC4020 plotter tape for microfilm production. It must 
be the last SC4020 subroutine call in the program in order to force 
emptying of any residual contents in the SC4020 output tape buffers at 
the end of the computer run. It is part of the SCORS subroutine package 
for the SC4020 plotter described in the publication "SC4020 Microfilm 
Recorder User's Manual" by Computer Sciences Corp. 

11. Subroutine EXIT 

A. Purpose — this is an IBM/360 Fortran Programming System subrou- 
tine. It is used to terminate program execution by returning control to 
the IBM/360 Operating System. See IBM/360 System Reference Manual 
"Fortran IV Language" (form number C28-6515-7) in Appendix C: 
Fortran- Supplied Subprograms. 

NOTE: All subroutines of the program which are not described above are used 

in generating the SC4020 plotter tape for microfilm production. They 
are part of the SCORS subroutine package for the SC4020 plotter de- 
scribed in the publication "SC 4020 Microfilm Recorder User's Manual" 
by Computer Sciences Corp. 

12. Program Input/ Output Details 

For a more in-depth explanation of input/output details than that following see 

Appendix J which contains a listing of the actual IBM/360 Job Control Language 
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used in the program. It also shows where the TI, TF selector cards go in the 
deck and their format. 

A. Inputs 

1) the data tape from the ATS-5 ground station 

a) read by fortran in the main program on unit 11 

b) tracks — 9 

c) mode — binary (as usual for 9 track tapes) 

d) density — 800 BPI (as usual for 9 track tapes) 

e) format — see Appendix B and section entitled "Ground Station 
Data Tape Format" above. 

NOTE: The RECFM (physical record format) parameter on the DD card for 

the input data tape was set to U (meaning undefined) because the original 
data tapes may have occasional short length physical records. 

2) the TI, TF selector program input cards 

a) read by fortran in the main program on unit 5 

b) position in deck — immediately after the GO.DATA5 DD * JCL 
card of the program deck in chronological order (see Appen- 
dix J) 

c) format — data fields other than "PRTSEL" are read by the 

I fortran field specification and must therefore be right-adjusted 
in the card columns assigned to the data field . 

data field main program card columns assigned to the 

fortran variable and description data field on the input card 

SYM, causes termination of computer 2, 3 

run if the characters XX are punched 
in the card columns assigned to it 

IYRI, last 2 digits of the year of the 
begin 
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time (TI) of the data interval selected 7, 8 

for processing by this input card, for 
example if the year is 1970, the char- 
acters 70 are punched 


IDYI, the day of year of the begin time 10, 11, 12 

(TI) of the selected data interval 

IHRI, the hour of the begin time (TI) of 14, 15 

the selected data interval 

MNI, the minute of the begin time (TI) 17, 18 

of the selected data interval 

ISECI, the second of the begin time (TI) 20, 21 

of the selected data interval 

IYRF, the year of the end time (TF) of 34, 35 

the selected data interval 

IDYF, the day of year of the end time 37, 38, 39 

(TF) of the selected data interval 

IHRF, the hour of the end time (TF) of 41, 42 

the selected data interval 

MNF, the minute of the end time (TF) 44, 45 

of the selected data interval 

ISECF, Ihe second of the end time (TF) 47, 48 

of the selected data interval 

PRTSEL, the variable on the input card 50, 51, 52 


used to select the time-ordered numeri- 
cal data value paper print-out display 
(see Appendix E for sample) of the data 
in the time interval selected by this TI, 
TF card. This display is outputted only 
if the characters PRT are punched in 
the columns assigned. If the assigned 
columns are left blank this display is 
not outputted. Variable PRTSEL is 
typed integer and read in by means of 
fortran field specification A. 
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NOTE: 


As stated above, when more than one TI, TF input card is used they 
should be put in chronological order as indicated in Appendix J. If they 
are not, the data selected by the TI, TF cards which are out of chrono- 
logical order with respect to the first TI, TF input card will not 
processed. However, the data selected by the first TI, TF card will be 
processed in any case if it is contained on the data tape mounted. 

NOTE: As explained in the section above entitled "Data Time Value Assign- 

ment" the TI, TF input cards must only specify data time intervals 
later in time than the ZERO YEAR value set in the program. 

NOTE: As mentioned above, to terminate the computer run properly, a card 

with the characters XX punched in columns 2 and 3 must follow the last 
TI, TF card. 

B. Outputs 

1) the tape which the Stromberg Carlson Plotter uses as its input to 
produce the microfilm consisting of the plot B display (see Ap- 
pendix F for sample) with the plot C display (see Appendix G for 
sample) . 

a) written by fortran in one or more of the subroutines of the new 
SCORS Stromberg Carlson 4020 plotter package on unit 10. 

b) tracks — 7 

c) mode — binary 

d) density — 556 BPI 

NOTE: In general several SC4020 plotter tape reels are outputted in a computer 

run since so many microfilm frames are necessary to display the data 
on a full ATS-5 ground station data tape. For such a full tape, i.e., one 
containing 30 days of data, approximately 7400 microfilm frames are 
produced from ten full SC4020 plotter tape reels outputted by the com- 
puter run. The new SCORS SC 40 20 plotter subroutine package and the 
IBM/360 Job Control Language DD card for fortran unit 10 shown in 
Appendix J in conjunction have all necessary logic for the program to 
request a blank tape to continue outputting SC4020 microfilm plotter tapes 
when a previous tape is filled up. 
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2) the optionally selected (see discussion of selector variable PRTSEL 
in the explanation of the TI, TF input card above in this section) 
time-ordered numerical data value paper printout display. 

a) written by fortran in subroutine ATSGPR on unit 6 

b) see sample of this display in Appendix E 

3) printed messages giving the beginning and end time (TI, TF) of 
each data time interval appearing on an input card read by the 
program; the beginning and end time of the mounted ATS-5 ground 
station data tape (the end time is printed only if the end-of-tape is 
encountered during processing of the data); the number of SC4020 
plotter microfilm frames residing on the plotter input tapes pro- 
duced by the computer run. 

a) written by fortran in the main program on unit 6. 
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APPENDIX A 


CONCEPT-LEVEL FLOW CHART 



FOLDOUT FRAME 


Read Next Input Card to Get Beginning 

and End Time (Tl, TF) cl New Date Time 2 1See No,e * 1 ond 2 B « low) 

Intervoi lo be Processed 


/D°et\ 

/ Card dull \ 

/ Read Hove the \ 
Signal (on XX) to End N 
the Computer Run? 


Print Beginning ond End Time (Tl, TF) 
of Data Time Interval Now Being 
Proceised an Output Listing 


NOTE: The term data block 
— means a physical record 

data bits on a tope 
bounded by speciol tape 
marks called “inter- 
record gaps". 



Indicate that the "End-ol- 

Run” Plots B <md C Are to COMMENT 7 

Be Dene 


NOTE: A "no" here 
tinuatioa al 


/Data BlociN. 
/Just Read Con-V 
train Data Within TimJy 
Interval Selected for 
Processing 5 
(Is TFREC Tl?) 


NOTE 1 : A number preceded by the word "commenl” appears 
next to certoin blocks in this (low chart. This is the 
number o( the comment in the progrnn listing in 
Appendix C which prefaces and locates the group 
ol progran statement! which actually accomplish 
the operations stated in the black. 

NOTE 2: Before reading this How chart, it is suggested thot 
the "plain language" general description of the logic 
How be read in the documentation under "Moin Pfogrom" 
on page 15. The comment numbers cited in this written 
explanation should indicate to the reader which parts 
ol the How chart correspond to the individual logicol 
functions ilemited on page 15. 


Get First Time (TIRED COMMENT W 
in Data Block Just Reod I — — — 


/ Does\ 

/Data BloeCV 
/ Just Reod \ 

' Contain Data WithinS. 
Time Interval Selected ' 
(or Processing 5 

(Is TIREC < TF?) 
















F0LD0U7 -FRAI^E 


Determine Best Ordinate 

Scale Limits for This 

Plot C from Max. and Min. 

of Averaged Data to be COMMENT 26 

Displayed in It (This Data ————— 





/ Signal to \ COMMENT 27 
r End Computer 
Run (an XX) Been N. 

Read from Input Card? 













Use Time of This Datum Value 
and Time Length Just Selected 
to Get End Time of Next Plot B 
Time interval Including Thi s 






APPENDIX B 


FORMAT OF THE ATS-5 GROUND STATION 
MAGNETOMETER DATA TAPE 





APPENDIX C 
PROGRAM LISTINGS 




2 FORMAT ( 20 0A4 .2 00 A 4 .200 A4 9 200 A4, 200 A4, 2C0A4, 200A4t 2 00A 4.200A4 ) 0000 0063 

4FH OAT 1 1 *"3 y^O VT 95“. OR i“T DA T C2Y37 .EQ -T9ST GITTO W3 WC00064 ' 

I RCI ND= I BLNK 00000065 

GO TO 704 OOCTOO OS'S 

C 00000067 

CZQM5 ***** ERR PRCCESSI NG T USE IDAT IF 2 DR MORE TIME FI ELDS WERE READ 0006S 
C * IF NOT CONTINUF READING UNTIL TRUE EVEN IF ERR FLAG ON - 00C00069 


To - ALT T5 ATS ““FRD BLOCK)' 


3 I R Cl NO - TA STRK" " " 

41 IFCIDATJ 1 i3) *NF • I 9S» AND • IDAT(2*3) »NE«I 9S) GO TO 704 


AU i 1 i 9 CINU=** / HlUAHJlUlJ=ltiOHl = <'llUU 

GO TO 41 00000075 

~C " "" “ 00000076 

CCOM6 ***** F ND OF TAPE PROCESSING Cl FILE PER TAPE) 00000077 

“C 000 00078 

4 CALL MSCLDPC TFREC. IYR2* IDY, IHR.MN.SEC) 00000079 

uf R I THTb * 7 06 ) IYR2f IDY. 1HR»MN.5EC 000 00080“' 

706 FORMAT C///1 X t 66HENC OUNT ERED END OF THIS PB TAPE - uA ST FIE-D DATA00000081 

I 'TTME ON TAPET IS T2 »1 H7Y IT»TH7»1 2 • IW.T2 rTH71Wi37 00000082 

45 IPCISTRTB.EQoC ) GO TO 1252 00000083 

“C “ '00000084 

CCOM7 ***** COMPUTER RUN DONF - JUMP TO PROCESSING OF DATA 00000085 


00000070 

00000071 

00CC0072 

00000073 


C-l 




TM2 = TIREC 

Z 

CCQM 1 A- ***** SEARCH THROUGH LOGICAL RECORDS OF BLOCK UNTIu NEXT 
C * TIMES ( TM 1 « TM2 ) DIFFERING BY 10 SECONDS OR LESS 


00000161 
OOCOC16Z 
OOOOC 163 
OOOOCI6A 
OOOOC 165 

OOOOC 167 


C-2 



LOGICAL RF CORD (STATION CODE, R CHANNEL USE CODE, YEAR 
“rODEr~TN“ARRAY IO” 


DO 16 1= 2 , I C 


IOFFST=4 

DO 25— r=-t V 3 - 

CALL PICK ( ID( I ).IDAT( 2, ISB-1 ), 0. I0FFST.4) 


CC0M19 ***** 
C * 


PUT THE 10 SETS OF H,D,Z,R DATA COUNTS OF THIS LOGICAL 
RECORD IN ARRAY I DCT - SKIP THE 48 DUMMY BITS AFTER 1ST 
H.O.Z.R SET 


CALL PICK ( IDCT< 1,1) , IDATT2, IS8-1T76, 16, 12) 

CALL PICK ( IDCT ( 2 , 1 ) • I DAT ( 2, IS8-1 ) ,0,23,4) 

CALL PICK! IDCT < 2 , TITT'ITaTI' 37 TiTTTTTSl 
CALL PICK< IDCT (3, 1 ) * I DAT ( 3, I SB- 1 ) • C • 8 • 12) 

CALL PICK ( IDCT (4 , 1 ) , I OAT 13, I SB- 1 ) , 0, 20, 12) 

I WD- 5 

NOBTS=l2 

NQTUSD=16 

DO'lTT= 2,10 “ 

DO 12 J= 1.4 
IWC = c 
NPTS0B=C 

14 CALL PICK ( IDCT( J, I ), IDAT ( IWD, iSB'-l j , IWC,NBTUSD*N0BTS) 
NBTS0E= NETS OB ••-NOB TS 
NBTUSb=NE TUSP+ NOBTS ~ 

I p ( NST SC B ,FQ* 12) GO TO 13 


00000225 

— OOCDC226 

00C00227 

— 

00000229 

00000230 

0000C231 

00000232 “ 

00000233 

00000234 

00000235 

00000236 

00000237 

50000238 

00000239 

00000240 

00000241 

00 0 O C 2 42 

00000243 

00000244 

000 00245 

00000240 " 

0000C247 


C-3 



0131 


NOBTS= 12 - NBTSOB 

1 1 1 1 1 1 1 1 — 

0132 


IWD=I WD+1 

00000249 

0133 


NBTUSD=0 

000 00250 

0134 


I WC= 1 

00000251 

0135 


GO TO 14 

00000252 

0136 


13 I F ( NBTUSD »FQ • 3 2 ) GO TO 15 

00000253 

0137 


NOSTS = 32-NBTUSO 

00000254 

0138 


IF( NOBTS • GT ,12) NOBTS=12 

00000255 

0139 


GO TO 12 

00000256 

0140 


15 NOBTS = 12 

00000257 

0141 


IWD* IWD + 1 

000002"58 

0142 


NBTUSD* 0 

00000259 

0143 


12 CONTINUE 

000002(16 


C 


00000261 


CCOM20 ***** STORE THE 1 0 H*D,Z,R SETS AND THEIR TIMES AND ID 1 Nr O 

00000262 


C 

* IN THE PLOT 0 ARRAYS- IF AN H»D,Z,R SET DATA TIME 

00000263 


c 

* EXCEEDS PRESENT PLOT S END TIME, PROCESS THE B ARRAY 

00000264 


c 

* DATA (AVERAGING =OR PLOT C>. DO PLOT B • THEN GET END TIME 

00000265 


c 




c 

* PLOT B 

00000267 


c 


00000268 

0144 


1= 0 

00000269 

0145 


24 I F ( TTEO . 10) GO TO 5TT0 

Uooocz r o 

0146 


1=14-1 

00000271 


c 


OOC 00272 


CC0M21 ***** AFTER SELECTING 1ST OF THE H,D,Z,R SETS *0R STORAGE IN 

00000273 


C 

* B ARRAY (1 = 1; JUMP TO tJ ARRAY END 1 1 Mb LJMPUIA11UN , h 1 . 

00000274 


c 

* IF THE STARTING-OUT SWITCH IS ON (ISTRTB=0> 

00000275 


c 


OUU002 f 6 

0147 


I F ( I STRTB • NE • 0 ) GO TO 17 

00000277 

0148 


TFP8 = C.DO 

000 002 f 8 

0149 


I STRTB* 1 

00000279 

0150 




0151 


17 IF(TSPL( I ) ,LE, TFPB)GO TO 601 

00000281 


c 


0000C232 


cc 

0M22 ***** PLOT B ARRAY END TIME EXCEEDED - COMPUTE 30 SEC 

00000283 


t 

* AVEtRAGC DATA POINTS FROM DATA IN B ARRAY AND STORE 

0000 0284 


C 

* THEM IN C ARRAY - GET END TIME (TAVL) OF 1ST 30 SEC 

00000285 


" 0152 “ 

0153 


ac>3 j=o 

TAVL — T IPB4-30 00 C • 
7029 DO 7005 K= 1,3 
CPSUM { K ) •= 0. 


COMPONFNT SUMMATIONS AND COUNTS TD START ( ICT IS 00000287 

“XTOnT^OF DATA POI NTS rN”"BARRAY‘ TO BE" PROCESSED! 01300028 8 

00000289 
OOOOOZW" 
00000291 


0154 

0155 


0156 

0157 


0158 

0159 


7005 CPCT ( K ) *0 • 

7019 IF(J.EQ.ICT) GO TO 7017 
JsJH 

IF(TMBARY( J) .GT.TAVL) GO TO 7017 


_ 
CC0M23 


00000292 

00000293 

00000294 

00000295 

0W0TJ296 

00000297 

00000299 
90COC300 
00000301 
“DOT'D O30Z - 
00000303 


***** STILL WITHIN THIS 30 SEC TIME INTERVAL - ADD THIS B 

* Array DATA < SELECTED BY SUBSCRIPT J ) TO JO' SEC SUMS 

* {DATA GREATER THAN 1950 IN MAGNITUDE IS ERRONEOUS AND 

IGNORED) 


DO 701 6 k=i , 3 

IF(ABS(BT(K, J) I.GT.1950.) GO TO 7018 
CPSUM ( K ) *CPSUM ( K ) *BT i K , J ) 


0160 

0161 


1)0000304 
00000305 
0 0 0 0 C30 6 
00000307 


0162 
01 63 


CPCT ( K ) =C PCT ( K ) + 1 . 

7018 CONTTNUE 


0164 

0165 


GO TO 7019 


00 C 00308 
00000309 


_ c 

CC0M24 

~c 

c 

~C 


***** END TIME FOR THIS 30 SEC INTERVAL EXCEEDED - SUMS FOR 

■* THIS AVERAGE POTFIT " ARE COMPLETE - USE TU UkJTA IN H,0,Z 

* COMPONENT AVERAGES -FIRST SEE IF ENOUGH ( CPCT MN> =»0I NTS 

FTAVE BEEN USED FOR RELIABILITY Uh AVERAGE POINT 


00000310 
OOC 0031 1 
0DDD031 2 
C000C31 3 


5166 “ 

0167 

0168 " 

0169 

0170 “ 


00000314 
0000031 5 


751 ? 

IF< JSB.EQ.O) JSB=2 

TDF*D APS ( TMBARY( J)— TMBARY ( J5BJJ 


‘00000316 

00000317 


000 0031 B 
000 003 1 9 
*011000320 
00000321 


CPCTMN* 10000 ,/TDF 

IP iCPCTl l t « L I .CPCTMN.ANU.tHtniit.LI i CP C f MN • AND ,CPC 1 l 3 ),LT,XPCTMN) 


1 GO TO 7008 

T 

CCQM25 

_ c 


***** 


AFTER OBTAINING A 30 SECOND AVERAGE DATA POINT JUMP TO 
C ARRAY END TIME (IFPC) OJMHU T ATI O N , E T C » IF I HE 


STARTING-OUT SWITCH IS ON (ISTRTC*0) 


IFOSTRTC.NE.OJGO TO 18 


00000323 

00000324 

0000C325 

00X300326“ 

00000327 


C-4 


ISTRTC=1 


18 IF( C TAVL - 15000* ) «LE ® TFPC) GQ TO 7025 


CCOM26 ***** PLOT C ARRAY END TIME EXCEEDED - GET PLOT C VERTICAL 


C * 30 SEC AVERAGE DATA POINTS IN C ARRAY TO BE PROCESSED) 


XMIN= 1.E6 


XMAX = 

DO 7066 L^ltlCTC 


IF(BTAV(K»L) *E Q« 9999 • ) GO TO 7066 
IE(BTAVCK*L) .GTaXMAX)XMAX=BTAV<K»L) 

IFCBTAVC KiL) a LT . X M IN ) XM I N=BTAV CK*L) 

7066 CONTINUE 

DO 7013 L= 1 • 6 

IF( XMlN.GEaSCLMI (L) • AND.XMAX aJ_E a SCLM2C L) ) SO TO 7012 


7013 CONTINUE 
7012 I VSC ( K ) =1 

CALL ATSGPC ( TMAV • BTA V • ISBST A* ICT C» I VSC • SCLM1 « SCLM2 • TI PC ) 

I FRCT = I FRCT ♦ l 
C 

CCOM27 ***** SWITCH NSWXX * 1 IF REMAINDER OF DATA IN C ARRAY AT 


END OF COMPUTER RUN HAS BEEN PLOTTED - SWITCH NSW130=1 
C * IF PLOT C DONE BECAUSE C ARRAY FILLED— UP (POSSIBLE ONLY 

C * TO RFSET COMPUTATION FOR NEXT TFPC IN THIS CASE) 

__ 

IF ( NSWXX *EQa l) GO TO 1252 


IF(NSW130 aEOal ) TFPC=OaDO 
C 

C COM 2 8 ***** INITIALIZE FOR NEXT PLOT C ARRAY «*• GET ITS ENTTTI MEi TFPC 
C * USING TIME OF 1ST AVERAGE POINT TO BE STORED IN IT 


7024 ICTC=0 
7011 TIPC=TFPC 

TF PC = T FPC + 3600000. 

I F ( ( TAVL* 150 00 # ) a GE a TFPC ) GO TO 7011 

N g wl 30 _ 0 

ISBST As= I DBAR ( 1 *J) 

7025 IF(ICTC.LTal3 0) GO TO 9998 
NSW1 30=1 

GO TO 9994 
C 

CC0M29 ***** COMPUTE RESPECTIVE HtOaZ AVERAGE VALUES “FROM 30 SEC 
C * SUMS AND COUNTS AND STORE IN PLOT C IF IT IS NOT FUuL 


C * -IF AVERAGE UNREL 

C 

9998 ICTC=ICTC+l 

DO 7026 K-1,3 

BT A V ( K .TCTCT = 9999 . ' 

IF ( CPCT ( K ) • LT • CPCTMN) GO TO 7026 


7026 CONTINUE 


7008 IF(J.EG.ICT) GO TO 7030 


CC0M30 ***** GET END T l MF (TAVL) OF NEXT 30 SEC AVERAGING INTERVAL 


00000329 
000 00330 
00000331 
0 000 0 332 
0000C333 


0000C335 

00000336 

00000337 

00000338 

00000339 


00000340 

00000341 


00000342 

00000343 

0OOOC344 

0000C34S 

“00000346 

00000347 


00000348 
, 0000C349 

OOOOC350 

00000351 

00000352 

00000353 


00000354 

00000355 

00000356 

00000357 

“00000358 


00000361 

“00000362 

0000C363 

O00DC364 

0000C365 


uuuucjoo 
00000367 
"00000368 
00000369 
“00000370 
00000371 
000003T2 
00000373 
’ 00000374 
00000375 
"00000376 
00000377 


oocoours 

00000379 

00000380 

00000381 

'00000382 

00000383 


00000385 
"00000386 
00000367 
'000 00388 " 
00000389 


GET NEXT AVERAGE POINT*ETC# 


19 TAVL =T A VL +3000 0 a 


00000393 
00000394 “ 
00000395 


00000396 

00000397 



APPROPRIATE - NOW DO PLOT B - GET PLOT B VERTICAL 00000399 

SCAlX"LTmTTS ' “FOR ' h; o »Z RESPECTIVELY "F I RST -"C ALL' AT5GPPT“ 00000400 
IF SELECTED TO PRINT-OUT DATA IN THIS B ARRAY 00000401 

— 000 00402 - 

7030 DO 30 J— 1 » 3 00000403 


XMIN=1 a E 6 

XMAX=0 a 

DO 3 1 K=1 .ICT 

I F ( ABS ( BT (JtK) ) a GE ® 19 50 a ) GO TO 31 


00000405 

00000406 

00000407 


C-5 


00000409 


0219 



IF ( BT ( J, K ) ,LT« XM I N ) XM I N =BT(.J»KI 


00000409 

0220 


31 

CONTINUE 


000004 I 0 

0221 



DO 32 K=1 *6 


0000041 1 

0222 



IF( XMIN. SE.SCLMl ( K) . AND. XMAX sLt.SCLMitKJ i Gu TO JO 


^T5t)0O04 1 2 

0223 


32 

CONTINUE 


00000413 

0224 " 


TO 

I VSC< J)= K 


00000414 

0225 



CALL ATSGPB< TMBARY • BT • I DBAR « ICT, I VSC , SCLM 1 • SCLM2 • T I PB « 


00000415 




1 TFRLTH ( iHSC) , I FLS. ltFLGl 


000OO41 6 

0226 



IFRCT= I FRCT ♦ 1 


00000417 

0227 



IF(PRTSEL.NE.PRTBCO) GO fO 2001 


0000041 8 

0228 



CALL ATSGPR( TMBARY , BT » R « I DBAR, IFLG, ITFLG, ICTJ 


00000419 

0229 

2001 

IFISVM.NETXXBCDJ GO TO B04 


^00000420 


C 




00000421 


— 1 1 1 1 1 II ■ !■ Ill ■ 


It TIM 


c 


4 OF COMPUTER RUN FORCE LAST PLOT C AND TERMINATE AFTER 


00000423 


c 


* IT IS DONE INSTEAD "OF INI T I ACT ZING~FT)R NE X I FLOT 3 


OO0QU424 


c 




00000425 

OTTO 



NSwxX= 1 


00000426 

0231 



GO TO 9994 


00C 00427 


c 




OOUUU42H 


CCQM33 4*4** INITIALIZE FOR NEXT USE OF B ARRAY - USE DATA 


00000429 


C 


4 SAMPLE TIME INTERVAL EXISTING FOR LAST DATA POINT PUT 


O0OOO43O 


c 


* IN PREVIOUS B ARRAY TO SELECT TIME SCALE (SELECTED BY 


00000431 


c 


4 SUB 5 CHIP T IHSL) rUH NLX1 PLU1 B — Gh I fcNU 11 ML 1 ! r ro J 

00000432 


c 


4 USING TIME OF 1ST DATA POINT TO BE STORED IN IT - RESET 

THIS 0433 


c 


W C umhu i A [ i On rTFFB I u ZERO ) TF PLOT □ TTME LENGTH HaS 




c 


4 CHANGED OR LAST PLOT B WAS DONE DUE TO FILLED B ARRAY, I. 

E 

•00000435 


c 


4 — SWITCH NSW 730 SET — TO I — (POSSIBLE ONLY - IF DATA SAMPLE 


00000436 


c 


4 TIME INTERVAL CHANGES WITHIN DATA IN LAST B ARRAY) 


00000437 


c 





0232 


804 

I CT = 0 


00000439 

0TT33 



1MSLSV— ifl SL 


00000440 

0234 



DO 20 J s 1 .5 


00000441 

023b 



IF ( \ spLST 7UT7I INT(J)} GO TO 21 


000 00442 

0236 


20 

CONTINUE 


OOCOC443 

023/ 



I HSC = 5 


00000444 

0238 



GO TO 1255 


00000445 

CT239 


2! 

IHMC = J-l : 


0000044o 

0240 



IF ( I HSC ,E Q * 0 }.I HSC= 1 


00000447 

0241 


i 25b 



000004 4 8 

0242 


60 2 

TIPB= TFPB 


00000449 

0243 






0244 



IF(TSPL( I ) .GE.TFPB) GO TO 602 


00000451 

0245 



NSW730=0 


00000452 


c 




00000453 


CCOM34 4444* STORE T I ME *H»D*Z,R«ST AT I ON CODE.R CODE, DATA YEAR CODE, 


00000454 


C 


4 ERR FLAG, AND TIMING FLAG IN APPROPRIATE 3 ARRAY- rt.D.Z 


00000455 


c 




■ 1 II 1 1 1 1 M 


c 


4 ONLY IF A FIELD READING - GO BACK TO GET NEXT OF THE 


00000457 


c 


* 10 H.D.Z.R SAMPLES OF THIS LOGICAL RECORD TO PROCESS 


0060 0458 


c 




00000459 

0246 


601 

IF( ICT.LT *730) GO TO 9970 


00000460 

0247 



NSW730=1 


00000461 

0248 



GO TO 803 


00000462 

0249 

9970 

I CT= ICT + 1 


OOC 00463 

0250 



TM8ARY(ICT)= TSPLCI) 


00006464 

0251 



DQ 23 J=1 ,3 


00000465 

0252 


23 

BT ( J . I CT ) = FLOAT ( IDCT< J, 1 >>* .376408 - 2000. 


00000466 

0253 



R ( ICT) = IOCT ( 4, I.) 


00000467 

0254 



IF( ID(2) .EQ. 1 .or7Td(2) .Eo.2 .OR . ID < 2J .Ea. 3 ) RTTCT 1 = R (TC T) * . 976408 


“OOCT0046ff 




1-2000 • 


00000469 

0255 



I DBAR ( 1 , ICT)=ID( 1 ) 


“OUOOT5470 

0256 



ID5AR(2,I CT ) — I D( 2 ) 


00000471 

0257 



IDBAR ( 3 , I CT ) — ID( 3 ) 


000 00472 

0258 



IFLG < ICT)=IRCI ND 


00000473 

0259 



I TFLG ( ICT ) = I YmTND 


— 0000 0474 

0260 



GO TO 24 


00000475 

0261 



END 


000004/6 











BEING PROCESSED INTO THE EQUIVALENT IN MI LLI SECONDS- SINCE-00000482 


THE TIME DATA - SEE APPENDIX B FOR DATA TAPE FORMAT 


SUBROUTINE ATS GRT ( IMD I « I WD2 • IY R* TM ) 
DOUBLE PRECISION tM 
DIMENSION ITMELI9) 

I0FF ST=0 
DO 1 1=1.8 

CALL PICK ( I T ME L{ I ) , I*D1 ,0 . IQFFST.4) 
1 I OFFST= I OFF ST + 4 


ID Y= I TMEL (1 ) ★ 1 0 0 + I T MEL (21*10 +ITMELI3) 


MN = I TMEL< 6 ) * 1 0 I T MEL ( 7 ) 



IF< IDY.LT • I .0R.IDY.GT.366.0R.I HR .GT .24. OR .MN.GT.60 .OR .ISEC.GT. 60) 00000502 


CALL MSZRDPI I Y R . I DY • I HR , MN , I SEC. TM ) 





***** “SUBROUTINE AT S GPR— GENER AfcS P^I NT- DUt ON TH£ SYSTEM OUTPUT 0000C51 1 

* UNIT OF THE DATA PRESENTLY STORED IN THE B ARRAY AS A 00000SI2 

* FUNCTION OF TIME - SEE APPENDIX E~FoRTA SAffP|_E OF THIS 00000513 

* PRINT-OUT 00000514 

06c 005T5 

SUBROUTINE ATS GPR < TMB AR Y , BT , R , I DBAR , I FLG , I TFLG , I C T ) 00000516 

DIME NS ION BT <3*730 ) . Rl730 TTlDB A« ( 3. 7 . 30) • IFL G < 730 ) , 1 St L8 L C 2 , 4 I , 000 0 0517 

l I RLBL ( 3*6 ) • I TF LG ( 7 30 ) 00000518 

DOUBLE PRECISION TMBARY("730l ' ~ " 00C00519 

DATA I STL BL/4HLYNN ,4HLAKE,4HTHQM« 4HPSUN , 4H» I NN » 4H I PZo ,4HT HE » 00000520 


14HPAS / * I RLBL / 4MN0T .4HUSE0.4H , 4HH AX.4HIS . 4H .4*0 AX, 00000521 

24HIS »4H * 4HZ AX.4HIS ,4H ,4HPROT,4HON E , 4HXP .4H0THE, 00000522 

34HR EX , 4HP 7 " “ " 00000523 

DATA LNCTR/0/ 00C00524 

COMMON/D A TE/MNTH, IOYMtH 00000525 

00 1 i« 1 , ICT 00000526 

IF(LNCTR.NE.O) GO TO 2 U0Q'ftC52T 

XSBSTA=IDBAP{ 1 . I ) 00000528 

ISBRc IDB AR ( 2 • X ) + 1 ” 0000C529 

WRITE (6. 3) ISTLBLU . ISBSTA) , ISTLBL(2* ISBSTA) • I RLBL ( 1, ISBR) • COO 00530 

— I tRL3L(2 * ISBR H IRL&LC3, ISBR) • lftLBL( IT tSBR") _ ,TRLBH 2 , ISBR) , 000 00531 

2IRL3L<3, ISBR) OOCOC532 

3 FORMAT ( l H i/// 9X, 54HATS-E MFM CANADIAN DOMINION 08SERV6 0 0 0 0533 

iATORY AT *2A4, IX, 43HM AN I TOBA MAGNETIC FIELD MEASUREMENTS 00000534 

2 R=. 3A4//10X , 95HDATE DAY OF TIME H AX IS D AXIS 00000535 

3 2 AXIS TOTAL FIELD FG MINUS R/8X, 93HYR MON DAY 00000536 

— 4 YEAS ffR m 5FC T5XMTMAST rGAMMAD) (GAMMAS) (GAMMAS) 00 C 00537 

5 PROTON ,3A4//) 0000C538 


93HYR MON DAY 
( GAMMAS) 


5 PROTON 
LNCtR= T.NC 


ARE INCORRECT, 


E~TD~^999V " IF ~A NY" OF THE '3 COMPONENTS 
GREATER THAN 1950. GAMMAS 


00000536 
00C 00537 
0000C538 
DC 0“0 053^5 
00000540 
00C 0054 1 
00C0C542 
00000543 
00 C 00544 


1 SORT (BT( 1 «I ) ** 2 ♦ BT (2,1) **2 + BT(3,I)**2> 00000546 

CALL M S C L D P ( TMB’ATfYTT T VI 7R , TDY~,THR ,~MN, SECT " 0000C547 

WRITE (6, 4) I YR , MNTH * IDYMTH, I DY , I HR , MN , SEC , BT ( 1 , I ) , BT ( 2, I ) • BT ( 3 , I ) ,0000054 8 

— rBTPrn > .tflgxtt.ttfug rn ~ 00000549 

4 ^OPMAT(8X ,12 , I X, A3, IX. 12, 3X, I3,4X, 12, 2X, 12.2X.F4, l , 3X , F 7 . 1 , 3 X , F 7 . I 000 00550 
1 ,3X,F7.T;'4X,F6 .1 7TZXTFTTT ,1x,ai, ix,ai ) OOudOSST - 

LNCTR= LNCTR+ 1 00000552 

RETURN " - 00C0C554 

END" 00000555 


C-8 


(INDIVIDUALLY) OVER THE CHRONOLOGICALLY NEXT DATA TIME 


00000561 


OF THE PLOT FOR A COMPONENT IS SELECTED FROM SEVERAL 


ACCORDANCE WITH THE RANGE OF DATA DISPLAYED IN THE PLOT - 00000565 

THE HORIZONTAL (TIME) SCALE IS C ROSEN FROM 5E VER7SL "05000565 

POSSIBLE SCALES FOR THE BEST DATA DISPLAY RESOLUTION IN 00000567 

ACCORDANCE WITH ThE DATA SAMPLING TIME INTERv AL'“F3R“THE OO 000568 

FIRST DATA VALUE STORED IN THE PRESENT B ARRAY - THE 00000569 




* HORIZONTAL (TIME) SCALE IS DONE IN THE MAIN PROGRAM - SEE 00000571 

* APPEND IX' F" FOR - '® -- S7^MPUE" _ PLITT B "" 00000572 

00000573 

— SUBROUTINE ATSGPB (TMBARY . BT, IDBAR, ICT, IVSC.SCLMl , SCLM Z.TI PB , — 00000574“ 

1 TFRLTH * IFLG • I TFLG ) 00000575 


1 SCLM 2 ( 6 ) » ISTLBL( 2*4) *MB( 3) ,MT(3) • I TFLG (730) 

DOUBLE PTjFCTSITjn tmb" yj T30 jtt TP 97TSCL 

DATA I STLBL/4HLYNN • 4HLAKE • 4HTH0M > 4HPS0N, 4HWINN.4HI PEG *4HTHE . 

j 4HpA 5 y 

OATA NST ART/C / * MB/672 *355 #38/ • MT/38 , 355 » 672/ » I BLNK/ IH / 
COMMOn/D ATE/mNtH , IDyMTH 
IF (NSTART.EO. 1 ) GO TO 100 

CALC I DF R MV ( * H TJVOTLL IB^ , *645* • r 2T.^ ,^2:279"**) : 

CALL C AMR A V ( 35 ) 

NSTAftT=1 - 

100 CALL FRAMEV 




00000577 
00 0 00578 
00000579 
00500550 
00000581 


***** CHECK FOR ILLEGAL STATION CODE BEFORE PRINTING 
LABEL TIN THIS “FRAME 


I F ( 1 DB A P < l , ITVGT ,'4T " GO T 0 105 
IS9STA = IDBAR (1,1) 


AT ION 00000589 

00000590 
OOC 00591 

00000592 ‘ 

00000593 


--KTWWW PX7TNT-ALX -LABELLING FOR THIS PLOT B 


10b CALL PRTNTVl 8, 8HMANITDBA, 470,10101 
CALL PR I NT V ( 6 • 6HPLQT B. 830,1010) 


00000595 

00000596 

00000597 

00000595 

00000599 


CALL PF I NTV ( 6 ,6HDAY = ,792,993) 

CA LL LTABLVT FLO AT ( I OY ) ,840*993, 3,1,3) 

CALL PR I NT V ( 10, 1 OH DATE = ,864,993) 

CAL L “TTABL.VT FLO AT < I DYMTH ) , 944 , 9 93 , 2-, 1 . 2 ) 

CALL PPINTV(3, MNT H ,960 ,993) 

C AlL l a b l v ( r lO A T TT7RTTT050 i"993T2"n ,'2 7 ' ' 

CALL PR I NT V ( 2 , 2HUT ,450,13) 

“CALL' PRTNTVt "5 .5HATS-E , 3 0 0 , 131 
CALL PR I NT V ( l 3, 1 3HDGM INI ON OBSY.850,13) 

CALL PFTNTVtl , 1HH, 27,788) 

CALL PR I NTV ( 1 , IHD, 27,472) 

CALL PPINTV( 1, 1HZ.27, 156) 

CALL APRNTV(C,— (2,6,6HGAMMAS,12,442) 

* **** GENERATE GRID AND PLOT FOR THE H DATA, THEN DO 

* THE D AND Z COMPONENT ON THE SAME MICROFILM = 

* COMPO NENT PLOTTED IS SELECTED BY DO INDEX J 

DO 5 J = 1,3 
I SUB= I VSC ( J ) 

CALL SFTMIV(42,8,MB( J) ,MT ( J) ) 

CALL GRID 1 V( 2 , 0 • , f FRLTH , SCLM 1 ( ISUB) ,SCLM2( ISUB) .TFRlT 
1 ( SCL M2 ( I SUB) -SCLM 1 ( ISUB) >/6.,0, 0,0, 1.0,4) 


***** DRAW TIC MARKS FOR HORIZONTAL (TIME) AXIS FOR THIS GRID 


DO 1 1=1 ,59 

IF(MOD( I V 10) .FQ.C) GO TO 1 
XI = FLOAT ( I ) *T FRLT H/60 • 
i F ( M OOT r75) .E Q .O ) " GO To 2 

CALL L I NF V ( NX V( X I ) , NYV( SCLM 1 ( ISJB ) ) , NXV( X I ) , NYV ( SCLM 
GO TO 1 

2 CALL LINFV(NXV(XI ) ,NYV( SCLM1 ( I SUB ) ) ,NXV( XI ) ,NYV( SCL Ml 
“1“ CONTINUE 


00000601 

00000602 

00000603 

00000604 

00000605 

■ 0 0 C 0 0 60 6“ 

00000607 
00000608 
00000609 
000 0061 0 
OOC 00 6 1 1 
00 00 0612 
00C0061 3 
00000614 
THE SAME F OROOO 006 1 5 
RAME - THE 00000616 
00000617 
0000061 8 
00C00619 
000 0062 6T 
0000C621 
H/6. , 00000622 

00000623 
00000624 
THIS GRID 00C00625 
0000 6 6 26" 
000 OC 62 7 
00000628 
OOC 00629 
00000630 
l ( I SUB ) ) + 8 ) 0CC00631 

000 00632" 

( I SU8 ) ) + 1 6 ) 00000633 

OOC 0 0634 
00000635 


C-9 




C 

c 


***** DRAW nc MARKS FOR VERTICAL (GAMWA' UNIT) AXIS' FOR 1 

* GRID 

rn i s 

00000637 


c 






0OO0O633 

0043 



DO 3 1=1*59 




00000639 

0044 



IFtMOOtl. 10).EQ.0> GO TO 3 




00000640 

0045 



XI = FLOAT! I )*( SCLM2 ( ISUB)— SCLM 1 ( ISUB ) )/60 • + 

SCLM1 ( ISUB) 


00000641 

0046 



TF [ MOD [ I . 5) .EO .0 ) 60 TO ~4 




00000642 

0047 



CALL LINE VtNXVtO. » .NYVtXI > .NXVtO . > + 8,NYV( XI ) > 




00000643 

0048 



SO TO 3 




000 00644 

0049 


4 

CALL LINEVtNXV (0. I .NYVCXI I .NXV <0.)+16. NYVtXI ) ) 




00000645 

0050 


3 

CONTINUE 




00000646 

0051 



DO 5 1= 1 .ICT 




00000647 


c 






015000643 


c 


***** COMPUTE HORIZONTAL COORDINATE OF THIS 

DATA 

PO INT 


00000649 


r 


^ (SELECTED BY DO I.NDEX I) - PLOT T AND/OR r 

FL AG IF 


000 0065O 


C 


* APPROPRIATE POR THIS DATA POINT 




00000651 


c 






(JUUU0662 

0052 



T= TMBARY ( I ) - TIPB 




00000653 

SOS'S 



IF { J • Nf"* 3 .Oft • I FLG ( I ) »EQ • I BLNK > GO TO *0 




00000664 

0054 



CALL PRINTV! 1 . IFLG! I ) ,NXV(T) ,NYV(SCLMl (ISUB) >- 

20) 



00000655 

0055 


46 

1F( J s n£« 3 «OR.I TFLSn ) • EO • I BLNK ) GO TO 40 1 




00000656 

0056 



CALL PRI NTV ( 1 • ITFLG( I ) ,NX«(T) , NYV(SCLM1( ISUB > ) 

-29) 



00000657 

00 57 


“431“ 

IX=NXV ( Y ) 




0000066b 

0058 



ARG2=BT( J # I ) 




0 0000659 

0039 



I Y = NY V l ARG2 ) 




OOOOOooO 


c 






00000661 


c 


***** PLOT VALOF'OF SfcLtCThU FILLS COMPUNtN! 

On 

IIS GRID 


OOU (JU662 


c 






00000663 

00 60 


5“ 

CALL PLOTV( IX* IY .42) — 




00000664 


c 






00000665 


c 


' ***** PRI NY THE TI ME ! HUUK* M INOTL *UR SECunD) 

HEPREb ENTET 

d T EATTruJ 00006 bo 


c 


* VERTICAL GRID LINE ON THE PLOT AT BOTTOM OF PLOT BENEATH 

00000667 


c 


* THE LINE 




0000C668 


c 






00000669 




ARG= SCLM 1 ( ISUB) 




00000670 

0062 



IY8=NYV( ARG) 




00000671 

010 63 



TSCL= TIPB ' “ 




000006 't 2 

0064 



DO 6 Is 1 «7 




00000673 




CALL mSCLOPITSCL. IYR , IDY . I NR. MN, SEC) 




000006 f 4 

0066 



T= TSCL - TIPB 




00000675 

0067 



I X=N X V ( T ) 




00000676 

0068 



IF ( I .FG. I *OR. I •E0.7,0R,TFRLTH,NE .60000. ) GO TO 

7 



00000677 

0069 



ISEC= SEC 




00000678 

0070 



CALL LABLVC FLOAT ( I SE C ) * IX— 8* I Y 8— I 2*2* 1*2) 




00000679 

007 1 



GO TO 6 




ODD 00660 

0072 


7 

CALL LABL V( FLO AT ( I HR ) • I X— 16*1 Y3— 4 * 2« 1 • 2 ) 




00000681 

0073 



CALL LABLVtFLOATt MN ) , IX. I YB-4. Z. ITS' J 




0oooo6az 

0074 


6 

TSCL = TSCL ♦ TFRLTH/6. 




00000683 

00 75 



“RETUPTT 




00000684 

0076 



END 




00000685 


C-10 







C ***** SUBROUTINE ATSGPC GENERATES A PLOT OF THIRTY 000 00690" 

C * SECOND AVERAGE H, D, Z* COMPONENT VALUES 00000691 
"C * (INDIVIDUALLY! OVER THE CHRONOLOG I C AXLY NEXT DATA TIME 00000692 
C * SPAN# I *E • CONTENTS OF PRESENT C ARRAY - THE VERTICAL SCALE00000693 
C * OF THE PLOT FOR A COMPONENT IS SELECTED FROM SEVERAL 0000069* 
C * POSSIBLE SCALES FOR THE BEST DATA DISPLAY RESOLUTION IN 00000695 
C * ACCORDANCE WITH THE RANGE OF DATA DISPLAYED IN THs PLOT - 00000696 
C * THE HORIZONTAL (TIME) SCALE IS SET AT 1 HOUR IN LENGTH 00000697 
C * AND BEGINS AT THE EXACT HOUR IMMEDIATELY PRECEDING THE 015W0698 
C * TIME OF THE 1ST DATA VALUE STORED IN THE PRESENT C ARRAY -00000699 
C * CHOICE OF THE VERTICAL SCALE FOR EACH COMPONENT AND THE 00000700 
C * HORIZONTAL (TIME) SCALE IS DONE IN THE MAIN PROGRAM - SEE 00000701 
C * APPENDIX G FOR A SAMPLE PLOT C 00000702 
C 00000703 


" SUBROUTINE A'TS GPC< TmCARY . ST , ISBSTA, ICTC. I VSC . SCLMXi SCU12, 71 PC) ' 
DIMENSION BT (3*130)* I VSC ( 3 ) • SCLM1 ( 6 ) t 


00000705 


0001 

0002 


1 SCLM2 { 6 ) • ISTLBLt 2*4) *MBC 3) *MT<3) 

DOUBLE PRECISION TMC ARY ( 1 30 ) , T IPC , TSCL 


06000706 

00000707 


DATA 1 STLBL/4HLYNN • 4HLAKE • 4HTHOM »4HPS0N«4HWI NN * 4H I PEoT 4HT *fE • 00000708 

14HPAS / 00000709 

DATA Mfi/672 • 355 ♦ 38/* MT /36T355r6‘72/T TFRLTH/ 3600000 ./ QUO 0071 0 

CQMMON/D ATE/MNTH* IDYMTH 0000071 l 


0004 


0005 

0006 
00 07 


CALL FRA"MEV 


00000712 

00000713 


***** CHECK FOR ILLEGAL STATION CODE BEFORE PRINTING STATION 
* LABEL ON THIS FRAME 


“000 00714 
00000715 


00000717 


0008 


I F ( I SB ST A • GT • 4 ) GOTO 105 


0009 

C 


CALL 

PRINTV(8,ISTLBL(1. I S 8 STAJ .200.1010) 

— 06000718 

00000719 


C 


***** PRINT ALL LABELLING FOR THIS PLOn 

00000720 


C 




00000721 

00 10 


10 b 

CALL 

PR 1 NTV { 8 » 8 HMANITQBA > 470 ■ 1 010 ) 

00000722 

0011 



CALL 

PR INTV(6*6HPL0T C.830.1010) 

00000723 

TOTS 



calL 

PR iNTV ( 20 , 2CHI 30 SECOND AVERAGES } « 270 * 993 } 

OOOOO / 24 

00 13 



CALL 

M SC LDP ( T I PC * IYR * IDY * IHR.MN. SEC) 

00000725 

0014 



CALL 

PRlNTVt - 8 ~. 6 HBAY 792, 9935 

00000/ 26 

00 15 



CALL 

LABLV( FLOAT ( IDY) ,840,993, 3.1,3) 

00000727 

<56T5 



CALL 

PRINTVC 10.J0H DATE — ,864,993) 

OOOOO / 28 

0017 



CALL 

LA3LV( FLO AT ( IDYMTH) ,944,993, 2* 1,2) 

00000729 

5oTb 



CALL 

PRtNTV< 3. MNTH.960.993) 

— 00000730 

0019 



CALL 

LA 8 LV( FLOAT ( IYP) , 1000,993*2, 1,2) 

000 00731 

6020 



“CAU7 

PRINT VI 2, 2MUT.AS0.13J 

00000/32 

0021 



CALL 

PR 1 NTV ( 5 *5HATS— E ,300,13) 

00000733 

S?T 22 



“CALL - 

TPTNTVT 1 3. ISHDOMINIQN QBSr.850.18) 

001700734 

0023 



CALL 

PR I NTV ( 1 • 1HH.27, 788) 

00000735 




CALL 

PR I NTV ( 1 . IHD. 27,472 ) 

ot j oou 735 

0025 



CALL 

PR I NTV ( 1 * 1HZ.27, 156) 

00000737 

00 25 



CAlL 

A PH NT v(0«— 12,0, oHGAMMAS , 12,4421 

000 00 738 


C 




00000739 


t 


***** GENERATE GRID AND PLOT FOR THE H DATA* THEN □□ THE SA Mb h UHOOU 00 740 


c 


* 

THE D AND Z COMPONENT ON THE SAME MICROFILM - RAME - THE 

00000741 


C 


* 


■ 1 1 1 m i ■■ 


c 




00000743 

0027 



DO 5 

J = l .3 

06000/44 

0028 



I SUB 

=IVSC( J) 

00000745 

0029 



” CALL 

SETMIVt 42 ,8 , MBI J ) * MT ( J) ) 0WO0746 

0030 



CALL 

GRID1V(2,C. , TFRLTH, SCLMi ( ISUB) ,SCLM2 ( I SUB ) , TFRlT H/ 6 • , 

00000747 




1 (SCLM2( ISUB)-SCLMU ISUB) )/ 6 ., 0*0*0, 1 . 0 .41 

OOOOC743 


C 




00000749 


c 


***** DRAW TIC MARKS FOR HORIZONTAL (TIME) AXIS FOR THIS GRID 

06660756 


c 




00000751 

003 1 



DO l 

1= 1 ,59 

00000752 

0032 



IF ( MOD (I* 10)*EQ*0) GO TO l 

00000753 

0033 



XI = 

FLOAT ( I ) * TFRLTH/ 60 • 

60560754 

0034 



IF ( MOD (I*5)«FObO) GO TO 2 

00000755 

0035 



CALL 

L INF VI NX V ( XI > ,NYV( SCLM 1 ( ISUB) ) ,NXV( XI ) , NYV ( 5CLM 11 I S JB1 )+ 8 ) 

■ i ri i n i — 

0036 



GO TO 1 

00000757 

003 7 

. 

— ? 

“CALL" 

“LI NF VfNXVCXr ) , NYV (SCLMI TT5UB) ) , NXV( XTT*NYVT SCLWI USUBD+16) 

00000758“ 

0038 


1 

CONTINUE 

00000759 


c 




00000750 


c 


***** DRAW TIC MARKS p QR VERTICAL (GAMMA UNIT) AXIS FOR THIS 

00000761 


c 


* 

GRID 

OOO 00 7 62 


c 




00000763 

0ff39 




“I =1 # 59 

00000764 

0040 



IF (MOD (I*10)*EQ»0) GO TO 3 

0000C765 


C-ll 




0042 



IF ( MOD ( I *5) o EQ oO } GO TO 4 

00000767 





CALL LIKIEVINXV 10.) .NYVIXI.J .NXVtO . J+8,NVV< Xl’J J 

0000076B 

0044 



GO TO 3 

00000769 

50^3 



CALL LINE V(NXV (0 « ) »NY VCXT)'*NXV (0 • ) 4- I6»NYV( X Z 1 > 

000 00 7/0 

00 46 


3 

CONTINUE 

00000771 

00"47 

C 


DO 5 T= TTTCTC 

000007 72 
00000773 


c 


*** ** C OMPUT E HORI ZflNTAC COORDI NA T£ CTF I HI 5 DATA P(J INI 

00000774 


c 


* (SELECTED BY DO INDEX I) 

00000775 


c 



00000/76 

0048 



T=TMCARY< D-TIPC 

00000777 

0049 



IX=NXVt i ) 

00000778 

0050 



ARG2=BT( J,I ) 

00000779 

0TT5T 

c 


IY=NYV(AHG2) 

00000780 

00000781 


c 


— mFrmp PLCJT VAL'tlE'OF' SELECTED FIELD C0HP0NEN1 ON ITS GRT3 

000 00 7 82 


c 



00000783 

0052 

c 

— 5” 

CALL PLOT VCT'X » IY • 42) 

00000785 


c — 


***** PR I NT TTTE TTME ( HUUK AND M I NUTE > REPREb tN lE j 3"T EAC 

H 0 0 0 0 0 7 66 


c 


* VERTICAL GRID LINE ON THE PLOT AT BOTTOM OF PLOT BENEATH 

00000787 


— c 


— * rHE LINE 



c 



00000789 

00 53 



ARG= SCLM 1 ( IbUB ) ' 

000 00790 

0054 



IYB*NYV( APG) 

00000791 

0035 



“T'SCL- — TTFC 


0056 



DO 6 1—1 *7 

00000793 








***** SUBROUTINE MSZRDP-CONVERTS A OA'tA TIME GIVEN BY YEAR* DAY 000008 0 6 

* HOUR, MINUTE. AND SECOND TO ITS EQUIVALENT IN A SINGLE TIME 00000807 

ZERO 

* IN DOCUMENTATION TEXT I 00000809 


00000810 

SUBROUTINE MSZRDP ( YR , DAY , IHR ,MN , I SEC , TMSZ YR » 00000811 

DOUBLE PRECISION TMSZYR, YMS, YINC , DAYMS.MSDY ” WO00B12 

INTEGER YR .DAY , YRSV 00000813 

common/zroyr/i zyr oOOOOei 4 

DATA YRSV/0/ 00000815 


00000816 

***** BY-PASS CALCULATION OF MILLISECONDS FROM ZERO YEAR TO YEAR00000817 
a OF THIS DATA TIME IF IT IS SAME AS YEAR” OF THE DATA TTWE ' 80000816 
* IN THE LAST CALL TO THIS SUBROUTINE TO AVOID UNNECESSARY 00000819 


COMPUT AT ION 


“08000820 

00000821 



l YMS=YMS-fY INC 


MSDY=IHR*3600000+ MN* 60000 ♦ISEC*1000 









MONTH* AND DAY OF MONTH (THE LATTER TWO OUTPUTS ARE IN THE 00000850 


QMMON SECTION NAM 

00000852 


UBROUTINE MSCLDP(TM, YR , DAY , HR , M I N, 

DIMENSION MNTHLMI 12) *MTHBCD( 12) 00000854 


DATA MTHBCD/3H JAN *3HFEB* 3HMAR * 3HAPR *3HMAY*3HJUN • 3HJU*. » 3HAUG » 00000856 


YR, YRORG* 


DAY, DYRM*HR*HRM 


COMMON/ZROYR/Y RORG 


TMwRK=TM 





MILLISECOND EQUIVALENT OF A YEAR CONTAINED IN THE INPUT 00000866 


FOR LEAP YEAR 


DO 1 1=1*5 



I F C M OD (YR*4) • E Q ■ 0 ) YRMS = 31 622400000 .DO 
~MW 

TMWRK=TMWPK - YRMS 


TIME WHEN ALL EXACT YEAR EQUIVALENTS IN MS HAVE BEEN 


00000878 

00000879 

00000880 


DYRM=DMOD( TMwRK, 8640 00 00 .DO) 000 00882 

000-00883- 

***** COMPUTE HOUR OF DAY FROM NUMBER OF MS REMAINING IN INPUT 00000884 


REMOVED - SIMILARLY GET MINUTE AND SECOND 



HR = DYR M/36 CO 00 0 


MIN= HRM/60000 


SEC= float<minrm)/iooo. 


***** FIND MONTH OF INPUT DATA TIME - TO DO THIS FIND THE MONTH 00000894 

-* HAV I NG A BEG I N AND END DA T O F^ T H E YE AR W HI C H 3R AClCE T TH E 0 0 0 00 8 95 

* DAY OF THE YEAR OF THE INPUT TIME AS FOUND ABOVE 00000896 


IF ( DAY «. GT .MNTHLM ( 1 ) ) GO TO 3 

IDYMTH=DAY 

MNTH=MTHBCD( 1 ) 


RETURN 


***** ADD 1 TO END DAY OF 


00000899 

00000900 


00000902 
00 0 00903 


JANUARY IF THE YEAR OF THE INPUT DATA TIME IS A LEA? YEAR 00000904 


3 I NC= 0 


I - ( YRMS . E Q . 3 1 6 22 400000 »D0 ) INC=I 
DO 4 1=2 , 12 

IF (DAY.LE •(MNTHLMI I )+lNC) ) GO TO 5 
4 CONTINUE 

***** ERROR RETURN 


00000906 


00000907 

00000908 

CT000C909 

00000910 

^50C009TT 

00000912 



IDYMTH=5C 00000914 

MNTH=MTHFRR 00 0 0 0 9T5~ 

RETURN 00000916 

oooowrr- 

***** CALCULATE DAY OF MONTH BY SUBTRACTING THE BEGIN DAY OF THE000009I8 


* DAY OF THE YEAR OF THE INPUT DATA TIME - IF THE MONTH OF 00000920 

* THE DATA - T I ME IS “FEBRUARY ! TS^EGInDAY OF THE YETAR TS HOO 0092 

* T HE SAME WHETHER THE YEAR OF THE INPUT DATA TIME IS A 00000922 

* LEAP YEAR OR NOT 00000923“ 

0000C924 


( I »EQo 2 ) IN 

IDYMTH=DAY— MNTHLMI 1-1 >-INC 00000926 

MNTrt=MTHBCD< I 5 0000092T" 

RETURN 00000928 


C-14 
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APPENDIX D 


"AUTO-FLOW" PROGRAM FLOW CHART 


NOTE: THE AUTOFLOW INSTRUCTION BOOKLET OBTAINABLE 

FROM THE APPLIED DATA RESEARCH CORP., OFFICE 
AT GSFC CONTAINS AN EXPLANATION OF THE NUM- 
BERS WRITTEN ADJACENT TO THE AUTOFLOW CHART 
BOXES. THESE NUMBERS ARE GENERALLY AUTOFLOW 
PAGE, BOX, OR FORTRAN STATEMENT NUMBERS. 



10/ 01-/ 70 


AUTOFl.Ua CHART SET 


C.S.F.C. ATS-S CRD STA MAC DATA PROG 


CHART TITLE - MAIN PROGRAM 


COM1 MAIN 

PROGRAM OF THE ATS-5 
GROUND STATION DATA 
PROCESSING PROGRAM 


■HITE TO 0£V / 

a / 

VIA FORMAT / 


I/YH = 69 
IYR = 69 


CCM2 ••••• READ 
NEXT Ti.TF UR XX IF 
NO KCHE 


| PH T SV = PHTSfc- 


/ READ FHOM OcV , 

/ S. ✓ 

/ VIA FUHMAT / 

/ 7J71 / 

/ INTO THE LIST / 


• LIST s SYM. IYRI • 
« IDY I • IHMI, UNI b 

• I SeCt. 1YRF. 

• ICYF, IHRF. MM- , 

• ISECF. PRTSEL 




/ 


FALSE . J , 


| PMTS&L * PHTSV | 


. A.OA. 

... is 


D-l 



10/0 IV 70 


AUTOFLOa CHART SET - G.S.F.C. ATS-S GRD STA 


DATA PROG 


CHART TITLE - MAIN PROGRAM 


/ 2000 

ca* 07 — >{ 


MS2RDP 
< 1 YH1 . 10Y1 . 
1HRI.MN1. 
IStCI.IU 


MSZHOP 
< I YRF, I0YF, 
IHRF.MNF . 
ISECF.IF) 


MSCLDP 

(T1 .1YR1.10VI. I 
IhRl.NNl.SECI) I 


MSCLOP I 
(TF.IYHF.1DYF. | 
1HHF.NNF.SeCF) I 


/ *RITE TO OfcV 
/ t> 

/ VIA FORMAT 
✓ 401 / 

/ FRCM THE LIST / 


LIST a 1 YR l • 

‘ 10 Y 1 • 1HMI. MN1. 

' SoCI. IVRF. IGYF. 
IMNF. MNF, ScCF 


COM3 44444 
CLIERMINE IF PRESENT 
BLOCK IS IN NE» Tl.TF 
• OR REAL) SLOCK 

IF AT TAPt BEGIN 
POINT (ENTER SEARCH 
PUCE) 


■ IdOTP .EG. 


COM4 •*•*4 BEFORE 
Ft AO INC NEXT dLOCK 
INTO CATA ARRAY ICAT 

• PNESET IT tG 

ALL S*S (THIS CAUSES 
9*S TO BE AT 

* ENG GF IUAT 

ARRAY IF bLUCK IS 
SHORT TO SHU4 ENG CF 
COCO 

a CA T A > 


/ READ FROM OEV , 
/ II / 

VIA FORMAT / 


/ INTO THE. LIST , 


LIST * 

( ( 1 GAT( J» 1 > . J a 

l.iai.i - v.ioi) 


COMS aaaaa ERR 
PROCESSING (USE IOAT 
IF 2 OR MORE TINE 
FIELDS MERE READ 
a IF NOT 

CONTINUE READING 
UNTIL TRUE EVEN IF 
ERR FLAG ON - 
a ASSIGN A FLAG 

OF F TO ALL DATA FROM 
BLOCK) 


| IRCIND * IASTRK | 


♦ -03. 25a > 


1 OAT (1.3) , 
I9S .AN) 
I OAT (2.3) c 
I9S 


/ READ FROM DEV / 
/ II / 

/ VIA FORMAT / 

/ 2 / 

/ INTO THE LIST / 


' ENO OF DATA? 4-a 


COMS *AA*4 FIND END 
TINE OF BLOCK (TFREC) 
IGNORING ANY 9*S IN 

* OAT A ARRAY 

(IBAT) - LSTLR IS 
SUBSCRIPT OF LAST 

a LOGICAL 

RECORD OF BLOCK - 
SUBA ATSGRT RETURNS A 
ZERO 

a IF TIME DATA 

LS BAD 


IDAT(I.I) • M 
I9S .AND. 
IDAT (2*1) »M 
19 S 


1 OAT ( 1 .3) .EG. 

I9S .UK. 

1CA T ( 2 • 3) .EG. 
I9S 


LIST a 

((IDATO.II.J > 
1.16). I a 2.101) 


a a TRUE 

a SEE NOTE a-a 
a ABOVE a I 


05.02 > l 


t .EG. 2 


1 IRC I NO = IBLNK 


D-2 



10/01/70 


aut of i.ob chart set 


S.S.F.C. ATS— S SRO STA MAS DATA PROS 


CHART TITLE - MAIN PROGRAM 


CCHt ***** END OF 
TAPE PROCESSING (1 
FILE PER TAPE) 


03.10*— >» 

^ 4 J 01 

1 1 MSCLOP H 

. J (TFREC.IYR2. H 

0 1 IDT. IHn.MN. H 

t | SEC) M 


/ *UIIfc TG OEV / 
/ 6 / 

/ VIA FORMAT / 

/ 706 / 

/ FRCM THE LIST / 


1 NOTE 03 
• • * • • 
• LIST « ITR2. IDT. * 
• IHR • MK. SEC • 


02.00 > j 

*5 * 04 


• 1STRTB .EQ. 


COM 7 ***** COMPUTER 

RUN OOME - JUMP 70 
PROCESS INS OF CASA 
• REMAIN! MS IN 

PLOT B ANO C ARRAY. 
THEN TERMINATE 


« 8*06. 
... 803 


1252 | 08 

/ WRITE TO OEV / 
/ 6 / 

/ VIA FORMAT / 

/ 1253 / 

/ FROM THE LIST / 

1 NOTE 07 
***** 
• LIST ■ ISRCT • 

| 08 

H 

PLTNO H 

H 




•"exit ~* 
CALL EXIT 


D-3 
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10/01/70 


AUTOFLO* chart set 


G.S.F.C. ATS-5 GRO STA MAG OATA 


CHART TITLE - MAIN PROGRAM 


/ 902 

03.29— >j 


ATS6RT H 

( 1 OAT ( 1 . 1 1 a H 

I0ATI2. II.1YR. H 
TFREC1 H 


« • TRUE 

TFREC .£0. «-» 

• 0.00 * I 


i LSIUN c 1 


C0M10 ***** BLOCK 
CMOS AFTER TI OR TAPE 
IS AT BEGIN POINT - 
COMPUTE 

• DATA BLOCK 
BEGIN TIME (TIREC1 - 
IF TIREC LESS THAN TF 

• SEARCH MODE 
IS COMPLETED - IFSTLR 
IS SUBSCRIPT OF 1ST 
GOOD 

• LOGICAL 
RECORD OF BLOCK 


03.07*— > 
BYO 


* TFREC .GE. TI 


ATSGRT 

(IDAT(l.I). 

IDATI2.I1.IVI 

TIRECI 


C0M9 ***** BLOCK 
PRE-OATES Tl.TF - 
ENTER OR CONTINUE 
SEARCH MODE - 

* LAST LOGICAL 
RECORD OF LAST BLOCK 
PROCESSED MUST BE 

* DISCARDED DUE 
TO BREAK IN TIME 
CONTINUITY IIPRELR 
SWITCH _ 

* 'SET TO ZERO) 


TIREC .EQ. 


IFSTuR = J 


NSCLDP 
ITI RFC. IVR2, 
IOY. HP. MS, 

Sr Cl 


TIREC .GT. 


COK1 1 ***** ANOTHER 
BLOCK IN THIS Tl.TF 
IS FOUND - PROCESS IT 


IPR=_R ! 
ISGTO : 


'RITE TO 3£V 


/ VIA =3RMAT 
/ 7ji 

/ from THE LIST 


COM1? •••»« BLOCKS 
ON ORIGINA. TA»S TF AD 
BY COPY PROGRAM WITH 
AN 

* t/3 ERROR 
HAVE NON-EERO HE < 
CHARACTER III 
INSERTED IN 

* 35T H HE* 

DIGIT OF C3*I ED 9.3CK 
- ASSIGN » 'LAG T3 
DATA 

* IN SUCH 
BLOCKS 


IRCIN3 * IAST». 


CONI? ***** LAS r 

L 33/ CAL PEC 327 3' 
-AST 3L 3C I PROCESS 
IS IN 1ST 

• I 3 WORDS 3 
I DAT ARRAY - PRJC? 
IT ON-Y 1= THERE I 
NO 

• 3RE4K IH T 
CONTINUITY T 3 N?<r 
DATA 3L DC< 


I I S3 s I I 


SB = IFSTLR 
T M2 = TIREC 


COM14 ***** SEARCH 
THROUGH LOGICAL 
RECORDS 3F BLOCK 
UNTIL NEAT 

• ONE BOUNDED 
BY 

VALID. CHRONOLOGIC ALLY 
ORDERED TAPE 

• TIMES 
(TN1.TM2) DIFFERING 
BY |A SECONDS OR LESS 

• AND WITHIN 
Tl.TF IS FOUND 


25*--> I 


COMlf *••«• 33 ON T D 

NEAT SL3CK IN TI.T® 

IF AL. LDSICAl 
RECORDS 0 s 

* THIS 3LDCK 

EXCEPT LAST ONE ARF 
PROCESSES - SAVE THE 

* -AST LOGICA. 

RECORD (At 9ESINN1N5 
OF |HT ARPAYl FIR 

* ' PROCESSING 
■ HEN NEXT f»-3C< 1 S 
ACCESS* D IF T HFR r IS 
NO 

* riM* sap 

BET KEEN 3LDCKS 

( I PR EL R SET TO 1 ) 


3-SIN DO -DO 3 


D-4 



t P/N1/TN 


AUTOF.J* CHART SET - G.i.'.C. ATS-S 3*0 STA HAS DATA 9X35 


CHART TITLE - .MAIN PROGRAM 


cam 6 *»*•* find 
NEXT GOOD LOGICAL 
RECORQ (SEE CONSENT 
1A FOR 

• OEFIN1TIDN OF 
• • GOOD * • > 

9S«2B — ->* 


ATSGRT 

( 1 DAT ( 1 .ISBI . 
10AT(2,I S3} . 
IYR.TN2I 




ATSGRT 

( IOATO .ISBI* 
1DAT ( 2 • I S3 ) • 
IYR.TNI) 


TNI «EQ* 0.D9 ' 


*«*»• A 5339 
TAI.TN? RATA ■HA* SEEN 
F3JN3 - PROCESS THE 

* H.3.Z.A SETS 
O 1 3 AT A ( _ 331 CA. 
REC3R3I BO JN D6 3 if 
THIS 

* TNI.TH2 - 
F!N> DATA SAMPLE TINE 
INTERYA. (TSPLST > 

* - ASS I GY A 
*LAG 0= T T3 AL. 10 
H.D.Z.R SETS IF 

tsplst 

* aiFFERS *RON 
ENGINEERING SPEC 
VALJES 


TNI 1/10.00 


TRJE • 

« TSPLST .GT. 

* 100P0. 




• * TRJE 
• TSPLST .EG. *-♦ 

* TtNTCI) * | 


| ITH1NO » ILTAT | 


• 7 . 02 . 

... 999 7 
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AUT0FL3B C<ART SET - GiSiF.Ci UTS-3 SRD STA HAS DATA PR 3 3 


IA/oj/70 

CHART TITLE - MAIN PROGRAM 




/ 9996 / 

06.13— > 

) Ol 


| 1TMIND = IBLNK | 
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10/01/70 


AUTOFL 3 A CHART SET 


- 3.3. = .C. ATS-5 S» ST A MAG DATA PROG 


CHART TITLE - MAIN PROGRAM 












07. 36 > 


* * TRUE 
TSPL1II «LE. »-« 

* T*P3 • I 


(FALSE . IP 


• « TAVL 

1 » .1 

• T- PC 


COM2 2 ***** PLOT B 
ARRAY END TIME 
EXCEEDED - COMPUTE 3' 
SEC 

• AVERAGE DATA 
POINTS FROM DATA IN B 
ARRAY AND STORE 

• THEM IN C 
ARRAY - GET END TIME 
(TAVL) OF 1ST 3) SEC 

• INTERVAL -ROM 
B ARRAY BEGIN TIME 
AND ZERO 30 SEC 

• COMPONENT 
SUMMATIONS AND COUNTS 
TO START IICT Is 

• COUNT 0= DATA 
POINTS IN B ARRAY TO 
BE PROCESSED) 

PA.05P >1 


7005 | 09 

CPCT(X) * 0. I 


7019 • 


J « EQ • ICT *— 


COM2 • ***** END TIME 
= QR THIS 31 SEC 
INTERVA,. EXCEEDED - 
SJMS *0R 

* THIS AVERAGE 
POINT ARE COMPLETE - 
USE TO OBTAIN H. D.T 

* COMPONENT 
AVERAGES -=IRST SEE 
IF ENOUGH I CPCTNN) 
POINTS 

* HA/E BEEN 
USED =OR RE.IABl.ITT 
DF AVERAGE POINT 

>3a.ii* — >• 

7017 | 15 

| JSB - J “ 1 | 

r 

• 13 

• • 

• • 

false * • 

_« JSB .EG. 0 • 



TRUE 

J JSB 

.1.1 

* 

19 

ry- = 

DAT3S C TMSARYI J ) - 
TMDARVI JSB)) 

CPCTNM » 
IP^OD./TDF 


NOTE 1 

CPCTI 1 
CPCTMN 
C*CT ( 2 
CPCTMN 
C»CT (3 
CPC 

• LT. 

• AND. 
•LT* 

• AND. 

• LT. 

MM ^ 


2' 


COMPS *•»•• PLOT C 
ARRAY END TIME 
EXCEEDED - GET P_3T C 
VERTICAL 

* SCA-E LIMITS 
FOR H.O.E 

RESPECTIVELY I ICTC IS 
COUNT OB 

* 30 SEC 
AVERAGE DATA POINTS 
IN C ARRAY TO BE 
PROCESSED) 

* -do plot c 


TRUE * 

. BTAVK.L) 

*. SO. ORR9. 


| XMAX = BTAVIK._) | 


TRUE * * 

* ABSI BT ( K * J ) ) 

*.GT. 1 9S0 . • 


CPSUMI K | * 
CPSUM(K) ♦ 
BT(K, J) 


COM2 3 ***** STILL 
■ITHIN THIS 30 SEC 
TIME INTERVAL - ADO 
THIS B 

* ARRAY DATA 
(SELECTED 0Y 
SUBSCRIPT J) TO 39 
SEC SUMS 

• (DATA GREATER 
THAN 1950 IN 
MAGNITUDE IS 
ERRONEOUS AND 

« IGNORED) 


COM2 S •••*• AFTER 
OBTAINING A 30 SECOND . 
AVERAGE DATA POINT 
JUMP TO 

• C ARRAY END 
TIME CTFPC) 
CDMPUTATIDN.ETC. IF 
THE 

• ST ART I N3- OUT 
S* ITCH IS 3N 
IISTRTC*0> 


FALSE * • 

< * BT AVI X 1 _ ) 1 

• .LT. X MI N • 


XMIM * 3TAVIK..I 


Pfl.27* > j 

7666 • 


TFPC x O.DD 
ISTRTC * 1 


I 9. Ill 
• •• 702 A 
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10/01*70 


AUTQFL3* C4ART SET 


G.$.?.C. ATS-S 3RD STA HAS DATA PROS 


CHART TITLE - MAIN PROGRAM 


00.33* >*<— 


XM IN a-GE. 
SCLMKL) .ANI 
XMAX .LE. 
SCLHZ(L) 


ATSGPC 
(TMAV.BTAV, 
ISBSTA. ICTC. 
IVSC.SCLM1 . 
SCLM2.TIPC) 


1 IFRCT a IFRCT ♦ 1 [ 


C0H27 ***** SWITCH 
NSWXX * l IF 
REMAINDER OF DATA IN 
C ARRAY AT 

* ENO OF 

COMPUTER RON HAS BEEN 
PLOTTEO - SWITCH 

NS* 1 30x1 

* IF PLOT C 

DONE BECAUSE C ARRAY 
FILLED-UP (POSSIBLE 
ONLY 

* IF DATA 

SAMPLE TIME INTERVAL 
CHANGES - TFPC SET TO 
ZERO 

* TO RESET 

COMPUTATION FOR NEXT 
TFPC IN THIS CASE) 


• NSWXX aEOa 


C 9.09* >* 

COM28 ***** 

INITIALIZE =0R NEXT 
PLOT C ARRAY - GET 
ITS END TINE (TFPC) 

* USING TTNE OF 
1ST AVERAGE POINT TO 
BE STOREO IN IT 

* -ISBSTA IS 
STATION CODE 
SUBSCRIPT FOR PLOT C 
LABELLING 


TIPC = T e PC 


TRUE * ITAV_ - i 

* 15000a) aSEa 

• TFPC ' 


P8a 23— 
7025 


— > 

15 


COM2 9 •*••• COMPOTE 
RESPECTIVE rl.Da Z 
AVERAGE VA.UES *R3N 
33 SEC 

* SOWS AN3 
COUNTS AN3 STORE IM 
PL3T C IF IT IS N3T 
FULL 

* -I* AVERAGE 
UNRELIABLE SET IT T3 
9999a 


| NOTE I a 

• BEGIN DO L33» * 

* 7326 < * la 3 * 



j BTAVIKa ICTC ) * | 

t 9999a 1 


C0M3I ***** DATA O" 
THIS 3 ARRAY AVERAGED 
- PLOT C WAS DONE I" 

* APPROPRIATE - 
NOW 30 PLOT 3 - GST 
PLOT 9 VERTICAL 

* SCALE LIMITS 
FOR N.D.Z 

RESPECTIVELY FIRST - 
CALL ATSGPR 

* !' SELECTED 
TO PRIMT-OUT OATA IN 
THIS B ARRAY 


7530 | NDTS 2B 

• *•••••**•• 

• BEGIN 00 LOOP * 

• 303*1.3 * 


XMIN * 1.E6 
XMAX * 0. 


t NOTE SO 

• •**«**•*•• 

• BEGIN DO LOOP • 

• si < * t. icr • 


~i 


TRUE * 

— — * CPCTIO aLTa 
• CPCTHN i 



FALSE 

21 

T BT AV ( K. ICTC I ■ I 

| CPSUMtO /CPCTIO | 

7026 

22 


* * N3 

END 3* DO * 

• L33»T * 


08.23- — > I 
7 30 8 * 


COM3 5 ***** GET SM> 
TIME ITAV.) OF NEXT 
30 SEC AVERAGING 
INTERVAL 

• FROM TINE 0 * 
1ST B ARRAY POINT TO 
BE I NCLUDE3 IN 
AVERASE- 

* GET NEXT 
AVERAGE P3I NT .ETC. 


• • TRUE 

• A3SI3T(JaKI) •- 

••GEa 1953. • 


• • FALSE 
* 3T1J.KI aGTa *-♦ 

• XNAX * I 


| XNAX « BTIJ.K) | 
* — — -• 

09.32 >1 


* * FALSE , 

3MJ.K) .UT. * >1 

• XMIN * 


t XMIN * BTIJ.K) | 


P9. 31* — >[< ♦ 


• LOOPT * 

* • 

• • 


TFPC = C.OC 


I ^swisn * l | 
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10 / 0 1/70 


AUTDFL38 CHART SET 


SaitF.Ci US-5 STD STA WAS DATA P«3S 


CHART TITLE - MAIN PROGRAM 


XMIN .GE. 
SCLMKKI .AND, 
XMAX .LE. 
SCLM2IK) 



YES 

30 3 

“J 


• * 

Ives 


06 


ATSGPB M 

(TMBARV.BT, M 
IDBAR.ICT. H 
IVSC.SCLMl. H 
SCLM2.TIPB, H 
TFRLTMtlMSC). H 
IFLG.ITFLG) H 


1 XF ~Z_l. 1 ' rRCT * L* 




ATSGPR 

I TMBARV • BT • R « 
ID BAR, IFLG. 
ITFLG* ICT ) 


COM3 3 ••••• 

INITIALIZE “OR NEXT 
USE OF a ARRAY - USE 
DATA 

* SAMPLE TIME 
INTERVAL EXISTING FOR 
LAST DATA POINT PUT 

• IN PREVIOUS 8 

ARRAY TO SELECT TIME 
SCALE < SELECTED BY 

* SUBSCRIPT 
IHSC) FOR NEXT PLOT 8 
- GET ITS eNO TIME 
ITFPB) 

• USING TIME OF 
1ST OATA POINT TO BE 
STORED IN IT - RESET 
TH 

* COMPUTATION 
ITFPB TO ZERO) IF 
PLOT B TIME LENGTH 
HAS 

• CHANGED OR 
LAST PLOT B WAS DONE 
DUE TO FILLED B 
ARRAY, l.E. 

* SWITCH NSW730 
SET TO 1 (POSSIBLE 
ONLY IF DATA SAMPLE 

• TIME INTERVAL 
CHANGES WITHIN OATA 
IN LAST B ARRAY) 


ICT « 3 
ItfSCSV »“IHSC 


BEGIN DO LDOP * 

20 J * I, 5 • 


* * TRUE 

* TSPLST ,LT, * 

* TINT! J) * 


20 


I »HS= ; I J 

10 , 1 6 * — >| 

1255 * 21 


COM3 2 *•**• IF 
PROCESSING OATA 
REMAINING IN B AND C 
ARRAYS AT END 

* OF COMPUTER 
RUN FORCE LAST PLOT C 
AND TERMINATE AFTER 

* IT IS DONE 
INSTEAD OF 

INITIALIZING FOR NEXT 
PLOT B 


I NSBXK * 1 | 
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CHART TITLE - NON -PROCE DURAL STATEMENTS 


o. 

o 

a 

a 

< 

j- 

< 

Q 

< 


in 











V 



A 

ft 
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a 










H 


« 

CM 

A 




C\ 











X 


X 

ft 

V 


< 


(Y 

r 




*> 






ft 


ft 

A 
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ft 
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ft 

k- 

s 
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10/01/70 


AUTOFL3# CURT SET - G.5.F»C» ATS-S GRO ST A MAS DATA PROS 


CHART TITLE - SUBROUTINE MSCLOP ( TM . YR. DAY , NR < MIN .SEC » 


/ MSCLDP / 


**•»* A I NO YEAR OF 
DATA BY ADDING 1 TO 
ZERO YEAR FOR EVERY 

• MILLISECOND 
EOUI VALENT OF A YEAR 
CONTAINED IN THE 
INPUT 

• DATA TIME! IN 
MS SINCF ZERO YEAR J - 
MAKE PROPER 
ADJUSTMENT 

• FOR LEAP year 


□EG I N OO LOOP 


• A 3 3 1 T 3 3*0 
DAY DF YEAT DF A__ 
MONTHS DF THE YEAR 
AFTER 

• JANJARY I s 
THE YEAR 3= THE IN»JT 
DATA TIME IS A .£!> 
YEAR 


* **•* COMPUTE DAY 
OF YEAR FROM NUMBER 
OF MS REMAINING IN 
INPUT 

* TIME WHEN ALL 
EXACT YEAR 
EQUIVALENTS IN MS 
HAVE BEEN 

* REMOVED 


***** COMPUTE HOUR 
DF DAY FROM NUMBER OF 
MS REMAINING IN INPUT 
* TIME MHEN ALL 

EXACT DAY EQUIVALENTS 
IN MS HAVE BEEN 


* REN3VSD - 

SIMILARLY GET MINUTE 
AND SECOND 


HR * 0YRM/3S - ''*A'»i* 

MODI OYRM*36-?HDi* | 
M IN = HRM/5330M 


***** FIN3 MONTH OF 
INPUT DATA TIME - TO 
30 THIS FIND THE 
MONTH 

* HAVING A 
BEGIN ANO HN3 DAY OF 
THE YEAR MM! CM 
BRACKET THE 

* OAY 0= THE 
YEAR OF THE INPJT 
TIME AS FQUNO A30VE 


NOT: 2* 


• CALC J.ATE DAY 
OF MDMTH AY 
SUBTRACTING THE 3E51N 
OAY □- THE 

• YEAR OF THE 
MONTH CONTAINING THE 
INPUT DATA TINE «RDM 
THE 

• 3 AY OF THE 
YEAR DF THE INPUT 
DATA TIM: - I* THE 
MONTH a e 

• THE DATA TIME 
IS C S3RUARY ITS 3E3IN 
DAY 3= THE YEAR 15 

• THF SAME 
MHETHER THE YEAR D B 
THE INPUT OATA TIME 
IS A 

• .:*> Y EAT DR 
NOT 


•••*• ERRDR RETURN 


JDYMTH = DAY 

NTH = MTHQCDI t I 


IOYMTH * 5? 
MNTH s MTHERR 
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io/ri/ 7 n 


AUTUFLDW C i A3 T SF I 


DS 


- G.i.-.C. ATS-S 5 >*> 5T* MAS DATA O, 


CHART TITLE - * ROUTINE TO RICK OUT CONTIGUOUS HITS — O N C • 


1*. I ».I »< I'M 
HALT I?.* 


ST03F At»C ADD'S 


STC t,S L S3»’ 

F|* SHlri L-fr 
SINGLE I NS TR t 


STC T.SlOTa 
FIX SHIFT LE F I 
OOJSLF IN 5TR, 


I C 

Check rn« DOU-l_e 
| wnwo-CONHiNi n; 


SL SR I 

3LL R. 
SHI-T LFF 


NEXT I 
ST 

STiIRfc OAT A AT 


CSTORF v r<S 

I • » 1 ?< 1 


END OF AUTOFLDN CHART SET 


0*1 |NPJT STATEMENTS POOCFSSS0 


FKrCJTIDN 


s;c 
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APPENDIX E 


SAMPLE OF THE NUMERICAL DATA VALUE 


PRINT-OUT DISPLAY OUTPUT 



ATS— E MFM CANADIAN DOMINION OBSERVATORY AT THOMPSON MANITOBA MAGNETIC FIELD MEASUREMENTS R=*NOT USED 



E-l 



APPENDIX F 

SAMPLE OF THE MICROFILM NON-AVERAGED 
DATA DISPLAY OUTPUT (PLOT B) 






APPENDIX G 

SAMPLE OF THE MICROFILM AVERAGED 
DATA DISPLAY OUTPUT (PLOT C) 





APPENDIX H 

DEFINITION OF THE H, D, Z DATA COORDINATE SYSTEM 



GEOMAGNETIC NORTH POLE 



H,D,Z Coordinates. On any spherical surface concentric with the earrn 
the H axis points to the geomagnetic north, the D axis to the geomagnetic 
east, and the Z axis to the center of the earth. 
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APPENDIX I 

IBM 1800 PROGRAM FOR COPYING 
ATS-5 GROUND STATION TAPES 

Written by Dave Fisher of IBM 
at Goddard Space Flight Center 



*APPENUIX I - IBM 1800 PRUGRAM PUR COPYING 

* AT S— b GRUUlML) ST AT I UN TAPtS - WRITTEN BY 

*UAVE PISHER UP IBM AT GUUOARU SPACE PLIGHT CENTER 

¥ 


BEGIN 

L OX 

L 1 

-3600 


L U 


HE X9 

LUUP 

STU 

L 1 

ARE A+360 1 


MUX 

1 

1 


MDX 


LUUP 


CAL l 


M AGT 


UC 


L 1ST 


L U 


LIST 


BSC 


Z 


MDX 


*-3 


LU 


L I ST + 6 


CMP 


PUUR 


MUX 


W R I 1 E 


MUX 


WRITE 


CALL 


M AGT 

*1 APE 

DC 


L 1ST M 


LU 


L 1ST m 


BSC 


Z 


MUX 


*-3 


CALL 


T YPEN 


UC 


L STYP 


CALL 


EXIT 

WRITE 

CMP 


SIX 


MUX 


WRT 


MUX 


WRT 


LU 


UNE 

*RECURU 




STU 


ARE A+9 

WRT 

call 


MAGT 


UC 


WL 1ST 


L U 


WL I ST 


BSC 


Z 


MUX 


— 3 


MUX 


BEGIN 

L ST YP 

UC 




UC 




BSS 


b 


UC 


/2U11 


UC 


MES 

SIX 

UC 


6 

Ui\lE 

UC 


1 

HE X9 

UC 


/Q9 t iH 


PILL BUPPER WITH HEX 9‘S 

READ TAPE RECURI) 

TE SI BUSY 

CHECK PUR EUP 

WRITE PILE MARK UN UUTPUT 
TEST BUSY 

I Y PE ENDING MESSAGE 

CHECK READ ERRUR 

YES-SET PL AG IN UUTPUT 

WRITE UUTPUT TAPE RECURI) 
TEST BUSY 

GU TU PRUCESS NEXT RECURI) 


1-1 



MtS DC MtS2“MtSl 

MtSl DMES 'KATS-t COPY JUb CUMPL fcTfc!) • t 

MES2 BES 0 

LIS l DC 

DC 

BSS 4 

DC 

DC /2000 

DC AREA 

L I ST PI DC 

DC 

BSS 4 

DC 

DC / 8U0 1 

DC 

WL I ST DC 
DC 

BSS 4 

DC 

DC / 400 1 

DC AREA 

BUUR DC 4 

AREA DC 3600 

BSS 3600 

END BEGIN 
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APPENDIX J 

LISTING SHOWING THE IBM 360 JOB CONTROL CARDS 
AND THE TI, TF CARDS FOR RUNNING THE 
ATS-5 GROUND STATION MAGNETOMETER 
DATA PROCESSING PROGRAM 



